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Beneficiation of Low Grade 


Saskatchewan Uranium Ores IV 


L. W. CRAWFORD’, BRAD GUNN’, S. 


A low grade pegmatitic uranium ore from Northern Saskatchewan 
has been ungraded by flotation by factors of from 3.8 to 5.1 with 85 to 
75% recovery of U.O,, which occurs in the form of uraninite. The 
collecting properties of a number of sodium salts of pure saturated and 
unsaturated fatty acids have been studied under a variety of conditions. 
The most successful beneficiation of the uranium ore was attained using 
sodium hexametaphosphate conditioner, sodium myristate collector and 
oleic acid frother. The addition of the conditioner and collector with 
the ball mill charge gave concentration ratios and recoveries comparable 
to those obtained with deslimed feeds and also resulted in a decrease 


in the consumption of frother. 


kEvious publications from this 

laboratory ‘1-2-3 4) have described 
attempts to beneficiate by flotation 
low grade pegmatitic uranium ores 
from the Charlebois Lake district of 
Northern Saskatchewan. The present 
work is a continuation of the flotation 
testing program and describes results 
which have been obtained when the 
sodium salts of a series of pure long 
chain fatty acids were used as collec- 
tors. The effect of various phosphate 
conditioners cn the flotation of ura- 
ninite has also been studied. 


Experimental 


A Fagergren laboratory flotation 
machine equipped with a 600 gm. 
capacity (overflow capacity 3600 ml. 
of water) heavy pyrex glass flotation 
chamber was used in all tests, The 
collection of froth was aided by hand 
skimming, with a one in. wooden 
spatula. Unless otherwise noted, the 
ball milling, desliming, conditioning, 
and flotation procedures were similar 
to those previously reported’), In 
most of the tests to be described the 
deslimed ore (95% minus 65 mesh ) 
Was conditioned in the flotation cell 


for five minutes, a specified amount of 
a phosphate salt was added as a dilute 
solution in water and conditioning 
carried on for another five minutes. 
The collector reagent was then added 
and after five minutes conditioning air 
was introduced to see if any flotation 
would occur at this stage. The only 
collectors that showed any appreci- 
able frothing under these ‘conditions 
were sodium oleate and large excesses 
of sodium myristate. With ‘the excep- 
tion of one commercial mixture the 
collectors used were sodium salts of 
pure fatty acids*. When pure acids 
were used as starting materials, a 
weight (usually equivalent to 2.7 Ib. 
of acid per ton of ore) was neutral- 
ized with the stoichiometric amount 
of sodium hydroxide in ethanol. The 
soap was added to the pulp as a sus- 
pension in ethanol. After collection of 
any froth which appeared, the air 
supply was cut off and measured 
amounts of a 5% solution of oleic 
acid in ethanol were added. After 
approximately two minutes condition- 
*The authors are indebted to B. M. Craig of 
the Prairie Regional Laboratory of the National 


Research Council for supplying the pure fatty 
acids. 
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ing, air was again introduced and the 
froth which formed was collected in 
a separate receptacle. More oleic 
frother was added and the additional 
froth again collected separately. This 
procedure was repeated until it was 
considered that sufficient concentrate 
had been collected to show the effect 
of the reagents used. 

In a few tests the desliming step 
was omitted or the conditioning and 
collecting reagents were added to the 
ball mill charge. In the latter case it 
was noted that even the longer chain 
fatty acid soaps (e.g. stearate) show ed 
a considerable tendency to froth with- 
out the addition of oleic acid. 


The ore sample, DE-53-1, used in 
this work has been previously de- 
scribed (3:5), It contained 0.05-0.08% 
of U,Ox, in the form of uraninite. 

All uranium analyses were by the 
fluorimetric method (5-5), Solid samples 
were heated for several hours with a 
1:1 mixture of concentrated sulfuric 
and nitric acids, diluted with water 
and neutralized by the addition of an 
excess of sodium carbonate ‘7)*. The 
mixture was then diluted to exactly 
500 ml., thoroughly mixed and al- 
lowed to stand for about two hours. 
An aliquot of the clear supernatant 
solution was then further diluted to a 
desired concentration with dilute nitric 
acid and 0.1 ml. aliquots of this solu- 
were taken for fluorimetric 
analyses. The total uranium recovered 
in the slimes, concentrate fractions, 
and tailings, normally agreed within 
+ 10% with that determined by ana- 
lysis of the feed. Material balances 
showed that 95-98% of the feed 
weight was recovered. 


tion 


©The effectiveness of this procedure for separ- 
ating uranium from carbonate insoluble con- 
taminants was checked by the authors using 
synthetic mixtures of uranyl and ferric salts. 
The method was found to yield reliable results 
within the accuracy of the fluorimetric analysis 
and within the concentration range encountered 
in this work. 
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Figure 1—Flotation of uraninite with soaps of saturated fatty acids. 


oo rr 
90+ 
80 + 
© 
oO 
© 
> 70+ 
iw 
oO 
60 + 
> 
a 
> 
50 + 
Oo 
oO 
40 + 4 
: / 
u Vi On 
O3s0+ / /// 
oc l[L/ 
lJ 
a 
20 + 
10 + 
oO ae nesiteill 
O 5 
Conditions: 
0.3 Ib./ton of ore of sodium hexametaphosphate 
3.0 lb./ton of ore of collector 
oleic acid frother as required. 
Legend 


A. Sedium myristate collector 
O  Sedium palmitate collector 
OQ Commercial soap collector 
©Q- Sedium laurate collector 
© Sodium stearate collector 


Regular pH determinations indi- 
cated that the pH of the feed pulp 
was 7.2-7.3. The pH of the pulp 
after the addition of collector varied 
from 8.5 to 9.5. The addition of un- 
neutralized oleic acid frother lowered 
the pH of the pulp slightly, the final 
value ranging from 8.0 to 8.5. 

The results of flotation tests are re- 
ported in graphical form (Figures 
-5). The percent of total U,Ox col- 
lected in a particular test has been 
plotted against the percent of the 
weight of the feed to the flotation 
cell that was collected. Since the 
slimes were usually removed prior to 
flotation the U,O, content and the 
feed weight used in calculating per- 
centages were the sum of the U;O, 
contents and the weights respectively 
in the individual concentrate fractions 
and tailings of each test. A point on 
the graphs therefore represents the 
cumulative percent recovery of the 


100 


U.O, fed to the flotation cell and the 
percent of the total feed weight col- 
lected in concentrates up to that 
point. The slope of any straight line 
drawn from the origin and intersect- 
ing a line showing the above relation- 
ship, represents the upgrading factor 
achieved. 

The slimes, when removed, usually 
contained 8-10% of the weight and 
20-25 of the U,O, content of the 
ore fed to the ball mill. It is believed 
that the amount of slimes could be 
reduced by grinding in a continuous 
ball mill-classifier circuit. 


Results and Discussion 


Figure 1 represents the results ob- 
tained when the sodium salts of lauric, 
myristic, palmitic, and stearic acids 
were used as collectors at the rate of 
3 lb./ton of ore on deslimed feeds. 
Sodium hexametaphosphate was used 


as a conditioner at the rate of 0.3 


20 25 


lb./ton of ore and oleic acid frother 
was added as required. Included 
this figure are the results of using a 
commercial soap* as a_ collector, 
Duplicate flotation tests showed good 
reproducibility. For the test condi- 
tions used, sodium myristate gave 
optimum recovery of the U;QOx values, 
e.g. 80% of the U,QO, recovered i 
15% of the feed weight. Sodium pal- 
mitate was almost equally effective 
as a collector for uraninite. Sodium 
laurate and sodium stearate gave high 
grade concentrates in the first frac- 
tions collected but did not continue 
to show as high a selectivity for 
uraninite as the myristate or palmi- 
tate did. The amount of oleic acid 
required to produce a suitable froth 
was observed to increase with the 
chain length of the collector mole- 
cule. For the test conditions noted 
the estimated amounts of oleic acid 
required to float 15% of the feed 
weights were as follows: 0.7 Ib./ton 
of ore with laurate; 1.1 lb./ton of ore 
with myristate; 1.4 lb./ton of ore with 
palmitate; and 2.5 lb./ton of ore with 
stearate. These results show that a 
sodium myristate rougher flotation of 
deslimed feed assaying 0.07% U;Ox 
vields a concentrate grading 0.56% 
U 30x, and containing approximately 
75% of the Us,Og in the feed in 10% 
of the feed weight. 

Figure 2 depicts the collecting pro- 
perties from deslimed feeds of the 
sodium salts of several other long 
chain acids. A typical test in which 
sodium myristate was used as collec- 
tor has been included for comparison. 
Only one point was obtained when 
sodium oleate (3.0 Ib./ton of ore) was 
used as collector as this reagent caused 
frothing without the addition of oleic 
acid. The other unsaturated soap, 
sodium erucate (3.0 lb./ton of ore) 
required only 0.3 Ib./ton of ore of 
oleic acid to float 15% of the feed 
weight. The unsaturated soaps al- 
though giving a total uranium re- 
covery comparable to that achieved 
with sodium myristate or palmitate, 
did not exhibit as high a degree of 
selectivity for uraninite. The sodium 
salt of azelaic acid (a nine carbon 
dicarboxylic acid) used at the rate of 
3.3 lb./ton of ore was definitely in- 
ferior to the straight chain saturated 
soaps as a collector for uraninite. The 
test with sodium acetate (1.5 lb./ton 
of ore of NaC,H;O,. 3H,O to give 
the same molarity of salt as in the 
standard myristate test) was done to 
show the collecting property of the 
oleic acid frother. To collect 15% of 
®*The commercial preparation consisted mainly 
of the sodium salts of palmitic and stearic acids 
with approximately 5% of sodium oleate. The 


authors are indebted to Dr. B. M. Craig for the 
analysis of this mixture. 
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the ‘ced weight 0.4 Ib./ton of ore of 
oleic acid was required. However, no 
selectivity for uraninite collection was 
exhibited by the oleic acid in the 
presence of sodium acetate. It was 
noted that the grade of the concen- 
trate tended to improve as additional 
amounts of oleic acid were added 
(see Figure 2). 

The results given in Figure 3 show 
the effects of variation in the type and 
amount of phosphate conditioner with 
myristate and stearate collectors on 
the flotation of deslimed feed. Three 
conditioners, sodium hexametaphos- 
phate (NaPO,),, sodium metaphos- 
phate (NaPO;),, and sodium dihy- 
drogen phosphate NaH.PO,, were 
used. The latter substance was tested 
because it is an hydrolysis product of 
sodium metaphosphate. With sodium 
myristate collector (3 lb./ton of ore), 
optimum results with respect to 
uranium recovery and concentrate 
grade were obtained with sodium 
hexametaphosphate used at the rate 
of 0.3 lb./ton of ore. The use of 0.45 
lb./ton of ore of this phosphate or no 
phosphate at all resulted in about the 
same upgrading factor which was 
lower than that obtained with 0.3 
lb./ton of ore of the phosphate. Sodi- 
um metaphosphate and sodium dihy- 
drogen phosphate appeared to have an 
adverse effect on the collection of 
uraninite with sodium myristate, the 
upgrading being lower than that ob- 
tained in the absence of phosphate 
conditioners. It is suggested that the 
sodium hexametaphosphate probably 
aids in the collection of uraninite 
through functioning as an effective 
dispersant. 

The use of sodium hexametaphos- 
phate (0.3 Ib./ton of ore) also im- 
proved the beneficiation of the ore 
when sodium stearate (3 lb./ton of 
ore) was used as the collector. 

The amount of oleic acid frother 
required to float 15% of the feed 
weight was found to be dependent 
on the type and amount of phosphate 
conditioner as well as on the collec- 
tor. This is illustrated by the data in 
Table 1. It is evident that the addi- 
tion of sodium hexametaphosphate 
conditioner results in a considerable 
increase in the consumption of oleic 
acid. This is in agreement with the 
suggestion that the phosphate acts as 
a dispersant thus increasing the 
amount of surface on which the col- 
lector or frother may be adsorbed. 
The sodium stearate being less soluble 
than the myristate might be expected 
to require more of the oleic acid to 
effect frothing as was observed. 

Figure 4 shows the effect on the 
flotation of uraninite from deslimed 


TABLE 1 
VARIATION OF OLEIC ACID REQUIREMENTS 


Oleic Acid to Float 
15% of the feed 
weight 


Collector Conditioner 


3 lb/ton of ore Substance Amt. Ib/ton of ore lb/ton of ore 


(NaPOs)¢6 0.45 1.0 


Myristate 

Myristate (NaPOs3)¢ 0.30 ol 

Myristate. “2 none 0.15 
Myristate (NaPOs3)n 0.45 0.30 
Myristate NaH.2PO, 0.45 0.14 
Stearate... (NaPOs3)¢ 0.30 ee 

Stearate... none Bo 
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Figure 2—Flotation of uraninite with soaps of saturated, unsaturated and dicar- 
boxylic fatty acids. 
Conditions: 
0.3 lb./ton of ore of sodium hexametaphosphate 
oleic acid frother as required. 


Legend: 
© Sodium myristate 3 Ib./ton of ore 
OQ Sodium oleate 3 = Ib./ton of ore 
© Sedium erucate 3 iIb./ton of ore 


A Sodium azelaate 3.3 lb./ton of ore 


_] Sodium acetate trihydrate 1.5 lb./ton of ore 
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feed of variation in the amount and 
method of adding sodium myristate 
collector. The feed was conditioned 
with 0.3 Ib./ton of ore of sodium 
hexamietaphosphate in all cases. Pre- 
liminary results had indicated that 
3 Ib./con of ore of myristate collector 
gave nearly optimum conditions with 
regard to both grade of concentrate 
and uranium recovery. When the 
amount of sodium myristate was in- 
creased to 6 lb./ton of ore and added 
as a single batch prior to aeration, 
16% of the feed weight was floated 
without the addition of oleic acid 
frother. It appeared that this amount 
of myristate was considerably in ex- 
cess of that required for optimum 
flotation of the uraninite and, although 
the final recovery of U;Ox was essen- 
tially the same as with 3 lb./ton of 
ore of myristate, the grade of the con- 
centrate was lower. When 3 lb./ton 
of ore of myristate collector were 
added in 1 lb./ton of ore increments 
(as indicated by arrows in Figure 4) 
together with sufficient oleic acid to 
cause frothing, both the grade of the 
concentrate and the final recovery of 
U,;0, were considerably lower than 
in the case where the same amount of 
collector was added as a single batch 
during the conditioning period. The 
amount of oleic acid required to float 
15% of the feed weight in this test 
was equivalent to 1.2 lb./ton of ore. 
(cf 1.1 lb./ton of ore of oleic acid 
when all of the myristate was added 
initially. ) 

Figure 5 compares the results of 
flotation tests in which the desliming 
step was omitted with those of stand- 
ard tests where the feed was deslimed. 
Sodium hexametaphosphate was used 
as a conditioner in all cases at the rate 
of 0.3 Ib./ton of ore. The sodium 
myristate or sodium stearate collector 
was used at the rate of 3 lb./ton of 
ore and oleic acid frother was added 
as required. Curve 1 represents the 
results obtained when sodium myris- 
tate was used in the standard proce- 
dure on a deslimed feed. Curve 2 re- 
presents the results obtained with the 
same amounts of phosphate and myri- 
state except that the reagents were 
added directly to the ball mill and 
the desliming step was omitted from 
the procedure. In this case the con- 
sumption of frother was considerably 
lower than that observed for the 
standard test (curve 1), 0.50 lb./ton 
of ore of oleic acid being required 
to float 15% of the feed weight as 
compared to 1.1 lb./ton of ore in the 
standard test. Curve 3 shows the effect 
of adding the phosphate and myristate 
reagents to an undeslimed feed in the 
flotation cell. About 1.1 lb./ton of 
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TABLE 2 
COMPARISON OF FLOTATION METHODS 


Method 1. Reagents added to deslimed feed 
























Reagent Ib/ton of ore 
Sodium hexametaphosphate See Paeeeay 0.3 
SOemsnne MY FIStAte. og cos sk oe ndk viecewes 3.0 
Oleic acid. ... sedhetiieg See ars 1.1 
_% of weight % U3;0Os oe recovery Ratio of 
fed to ball mill in fraction of U;0s enrichment 
slimes ; 7 0.169 23 3 
am oe lees | SM tea | Btls im 
rougher conc... 14.3 0.244 61.5 
tailings 77.9 0.011 bon 
Head (cale.). . 100.0 0.0567 100.0 
Method 2. Reagents except oleic acid added to ball mill 
Reagent Ib/ton of ore 
Sodium hexametaphosphate ............... 0.3 
Soca: HVRIBIATE. ©... 5 once ns bce cies ara 3.0 
CUS AGM res as lw do eas Caen 0.8 
_& of weight &% of U;0s ©] recovery Ratio of 
fed to ball mill in fraction of U;0s enrichment 
rougher conc... . 22.7 0.287 83.2 3.67 
tailings 77.3 0.017 16.8 
Head (calc.). . 100.0 0.078 100.0 
TABLE 3 
COMPARISON OF FLOTATION METHODS 
Method 1. Reagents added to deslimed feed 
Reagent Ib/ton of ore 
Sodium hexametaphosphate ............... 0.3 
Sodium myristate...... nivaeeanalaaes 3.0 
Oleic acid cae ee re 0.8 
“~ of weight % UOs Co recovery Ratio of 
fed to ball mill in fraction of U30s enrichment 
slimes. . 7.8 \ e 0.169 20? 23.4 \ we 5.17 
rougher conc.. . 6.7 14.5 0.441 0.292 51.6; />-9 
COU. occ seo: 85.5 0.0165 25.0 
Head (calc.). . 100.0 0.0567 100.0 
Method 2. Reagents except oleic acid added to ball mill 
Reagent Ib/ton of ore 
Sodium hexametaphosphate § ohceas Soe 0.3 
Sodium myristate... . oats san 3.0 
Oleic acid... . Pez ; 0.4 
/ of weight % U3Os () recovery Ratio of 
fed to ball mill in fraction of U;0s enrichment 
rougher conc.. . 14.6 0.400 75.0 SS Fy 
tadings.... 6. 85.4 0.0228 25.0 
Head (calc.). . 100.0 0.078 100.0 
ore of oleic acid was required to float added to an undeslimed flotation cell 
15% of the feed weight in this case. feed (curve 3). Curve 4 represents 
It is evident that much more efficient a test in which sodium stearate col- 
use is made of the reagents when they lector was used in a standard test on 
are added directly to the ball mill a deslimed feed, while curve 5 shows 
charge (curve 2) than when they are the results obtained when the same 
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Figure 5—Flotation of uraninite with sodium myristate and sodium stearate. 
The effect of adding the collector to the ball mill. 


Conditions: 


Conditioner: 0.3 lb./ton of ore of sodium hexametaphosphate 


Collector: 


Legend: 


Curve 1, < 
Curve 2, | 


3.0 Ib./ton of ore of myristate or stearate 
Oleic acid frother as required. 


Sodium myristate added to deslimed feed 
Sodium myristate added to ball mill 


Curve 3, O Sodium myristate added to undeslimed feed 


Curve 4, /\ 


Sodium stearate added to deslimed feed 


Curve 5, \/ Sodium stearate added to ball mill 


amounts of sodium hexametaphos- 
phate and sodium stearate were added 
to the ball mill charge and the deslim- 
ing step omitted. Again, ball milling 
of the reagents with the feed resulted 
in a considerable decrease in the 
amount of frother required to float 
15% of the feed weight, 0.5 Ib./ton 
of ore of oleic acid being used as 
compared to 2.5 Ib./ton of ore in the 
standard test on deslimed feed. Ball 
milling of the myristate or stearate 
with the charge evidently results in 
much more efficient dispersion and 
adsorption of the collector on the 
mineral particles so that in spite of 
the presence of slimes the total amount 
of reagents required to float a given 
fraction of the feed weight is reduced. 
When the stearate collector was 
added to the ball mill (curve 5) the 
grade of the initial flotation concen- 
trate was not as high as in the case 
where the collector was added to the 
deslimed feed (curve 4), but the final 
U;0, recovery and_ concentrate 


104. 


grade were higher and comparable 
to the results obtained with myristate 
collector. 

A more direct comparison between 
the results obtained by adding the 
reagents to the ball mill ‘charge and to 
a deslimed feed is shown in Table 2. 
Combining the slimes fraction with 
the rougher concentrate (Method 1) 
the total U;O, recovery would be 
84.9% in 22.1% of the head weight 
giving an upgrading factor of 3.84. 
Much the same result is obtained by 
adding the conditioner and collector 
to the ball mill (Method 2), the 
U;,O, recovery being 83.2% in 22.7% 
of the head weight, giving an upgrad- 
ing factor of 3.67. In addition to the 
elimination of the desliming step, 
Method 2 has the advantage that the 
amount of oleic acid required is lower 
and it is probable that the amount of 
collector could also be reduced be- 
low the figure given in Table 2 in 
actual practice in which a classified 
feed would be used. The presence of 


slimes usually increases the consump- 
tion of flotation reagents, but in this 
case the more effective dispersion of 
the collector during the ball milling 
has more than offset this increase. 

It would be possible to take off a 
higher grade rougher concentrate 
with the sacrifice of U,Og recovery, 
For example, assuming that only 75% 
of the U,Ox, values were taken in the 
slimes plus rougher concentrate 
(Method 1) or in the rougher con- 
centrate (Method 2), a ratio of en- 
richment somewhat greater than 5.1 
could be achieved. The figures shown 
in Table 3 indicate the results to be 
expected from the two methods for 
the above condition. In this case the 
tailings could be scavenged down to 
at least the grades listed in Table 2 
and such scavenger fractions recycled 
with fresh feed. Such a method 
should be tested on a pilot plant scale, 

Some preliminary attempts to pro- 
duce a high grade cleaner concentrate 
from the rougher concentrate frac- 
tions have not proven to be very 
successful as the whole of the rougher 
concentrate tends to float again even 
without reagent addition. This prob- 
lem together with a study of the ex- 
traction of uranium from _ flotation 
concentrates by leaching will be con- 
sidered in subsequent reports. 


Conclusion 

Low grade pegmatitic ore contain- 
ing uraninite can be upgraded by 
flotation by a factor of 3.8 to 5.1 with 
U;O, recoveries varying from 85 to 
75%. The most satisfactory combina- 
don of reagents involves sodium 
hexamet aphosphate as a conditioner, 
sodium myristate or palmitate as col- 
lector and oleic acid as a frother. The 
addition of the conditioner and col- 
lector directly to the ball mill charge 
decreases the amount of frother re- 
quired, 
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A Mathematical Study of Crystal Growth 


in a Cascade of Agitators 


J. N. ROBINSON? and J. E. ROBERTS® 


General equations are developed which describe, under certain 
specified conditions, the behaviour of crystals growing in a cascade of 
agitators. The crystal size distribution function for a system containing 
any number of agitators is derived and used to compare the screen size 
distributions of crystals grown in multiple and single agitator systems. 
It is shown that as the number of agitators in the system is increased the 
size distribution of the crystals becomes more uniform but the mean size 
decreases for a given rate of production per unit volume. 


Some applications of these equations to operational problems 


are shown. 


Introduction 


THEN crystalline materials are 
produced for use in the fertili- 
zer trade it is desirable that the size 
distribution be as uniform as possible. 
An undue amount of either very fine 
or very coarse material will make the 
product unacceptable to the con- 
sumer. Normally ammonium sulphate 
is produced in either saturators or 
vacuum crystallizers operated as single 
units and its size distribution, while 
acceptable, is by no means uniform. 
This study was undertaken in an at- 
tempt to determine whether an exist- 
ing cascade of agitators at The Con- 
solidated Mining and Smelting Com- 
pany of Canada Limited should be 
operated as ammonium sulphate sa- 
turators in series or parallel to pro- 
duce an acceptable product at maxi- 
mum possible rate. 


Previous mathematical work on 
crystal size distribution ‘':?) appears 
to have been confined to single vessel 
systems and this paper is an attempt to 
derive equations for crystal size dis- 
tribution in the general case of “K” 
agitators operated in series. The in- 
vestigation of retention times of dis- 


crete particles in a series of agitators 
is important in this problem and has 
a parallel in the determination of 
purging times of a series of con- 
tinuous flow stirred tank reactors. 
(3, 4,5), However, the fact that in a 
crystallizer we are dealing with a 
finite number of particles per unit 
time means that the time to reach 
steady state oper ration is finite whereas 
in the purging of a series of reactors 
the time for complete purging is 
theoretically infinite. 

The mathematical development falls 
into three sections; the description of 
particle retention time in a cascade of 
agitators, the equations describing 
crystal growth, and the combination 
of these two results to describe the 
final particle size distribution of the 
crystals grown in the system. 

The mathematical theory is shown 
to agree well with available operating 
data for ammonium sulphate. The 
treatment in this paper is in no way 
exhaustive, the authors’ intention 
being to suggest mathematical tools to 
assist further investigations. A num- 
ber of implications of this study and 
some areas for further study are in- 
dicated. 
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1. Development of a general expres- 
sion to show the retention time 
distribution of discrete particles 
passing through a system of “K” 
agitators of equal size. 


(a) Assumptions and Notation 


Let us choose a time scale with 
intervals which are discrete and 
arbitrarily small. Let the system con- 
sist of “K” agitators each of volume 
“V” and assume that a flow of “ 
per unit time enters and leaves each 
agitator in the system. Let “n” par- 
ticles per unit time be added to the 
first agitator of the system. The pro- 
blem at this stage is to determine the 
distribution of retention times 
amongst the “n” particles per unit 
time which leave the last agitator as- 
suming that an equilibrium state has 
been reached and that there is perfect 


mixing in each agitator. 

We will write u/V = r throughout, 
and it is obvious that 1/r is the mean 
retention time in any one agitator. 
The time scale may alw ays be chosen 
so that r<l. 


(b) Single Agitator Case 

Initially we will consider the reten- 
tion time distribution for a single 
agitator. 

Consider the “n” particles which 
enter the agitator in the time interval 
just prior to t = 0. In the first time 
interval, t?, the expected value of 
the fraction of these particles passing 
out of the tank is u/V =r. 

The expected number of initial 
particles leaving the tank in time in- 
terval t? is given by 

y? = nr 

In the interval t} : 

y! =r(n — or) = or(1 —r) 
and in general: 


ot = eee ~ ap 
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This series ends when there are “r”, the expression is just less than 


none of the original “n” particles left 
in the tank. Let us assume that t™! 

is the time period during which the 
last of the original particles passes out 


of the tank. We can therefore write: 


O<n nr {1 +(1—r) +(1-—r)? +4 
+ (1—r)™" |] <1 


which reduces to 


The required value of “m” is the 


least integer for which, given “n” and 
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for all 


aS 


one. Since practical cases 
r<1 and such a value must 
exist and hence “m” is determinable. 


One other bound which is readily 
determinable is the point in time after 
which less than one particle per unit 
time of the original “‘n” is expected to 
leave the system. Since we are dealing 
with discrete particles this means that 
after this particular time, intervals 
when one particle leaves the system 
will be interspersed with intervals 


when no particles leave. This bound 
is given by: 


md as 
Ya on 


nr (1—r)e'< 1 


The retention time distribution of 
the original “n” particles is, of course 
given by: 
nr (1 ry! (1) 
and it is obvious that this is the reten- 
tion time distribution of every “n” 
particles leaving the system per unit 
time after equilibrium is reached. 


y= 


When starting with a tank of fluid 
containing no particles equilibrium 
will be reached “m” time units after 
the addition of “‘n” particles per unit 
time begins. 


An example of the retention time 
distribution for a single agitator is 
shown in Figure 1. In this example 
values of n=100,000 and r=0.01 
were chosen, and it is found that 
m = 1148 time units and q = 688 time 
units. 


(ec) Multiple Agitator Case 


Here we will derive the retention 
time distribution for the case of “K” 
agitators of equal volume “V” oper- 
ated in series. 


“ 


We have shown that of “n” par- 
ticles under examination, a number 


yt-! = nr (1 — r)! are expected to be 
retained in a single agitator “t” units 
of time. 

We may also say that the prob- 
ability of one particle remaining in a 
single agitator “t’ units of time is 
given by: 


1 


Pr. (t) = r(1 —r) 

If we now consider a cascade of 
“K” agitators of equal volume “V” 
each with an inflow and outflow of 
“u” per unit time and a feed of “n” 
particles per unit time to the first 
agitator of the series, then under 
equilibrium conditions “n” particles 
will pass into and out of each agitator 
per unit time. Each agitator behaves 
exactly the same as the single agitator 
system. If we now consider one group 
of “n” particles which enter the first 
agitator of the cascade, we want to 
determine the distribution of reten- 
tion times among these “n’’ particles 
when the last particle leaves the cas- 
cade. Since each agitator behaves 
independently, we can write: 

The probability of a single particle 
residing in the first agitator a time 
“p,” Is: 


Pr (pi) = r (1 r),~! 


The probability that it will remain 
in the second agitator a time “p,” is: 


Pr (ps) = r (1 r)Py~! 
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id 


id the probability that it will re- 
main in the Kth agitator a time “px” 
iS: 

Pr (px) =r(1l —r)x- 

Now the probability of a particle 
remaining in the Ist agitator a time 
“p) is completely independent of the 
probability of it remaining in the 2nd 
agitator a time “p,” or of remaining 
in the Kth agitator a time “px”. 
Therefore the probability of the par- 
ticle remaining p, in the Ist, p, in the 
2nd, . . - px in the Kth agitator is 
given by: 

Pr (pi, pe px) =r(1 - r)! 


D1 : - \?K 


A(l — py" 

, (y+, 4 ‘a ~~. 

rk (1 — r) - 

However we are only interested in 

the distribution function for the total 

retention time “R” of the particle 

in the whole system where R = p, + 
DoF Bie 


But for any given R>K and for 
K>1 there is more than one way of 
assembling a total retention time of 
“R” from components p;, Pp. --- Px: 
since it does not matter how much of 
the retention time is accumulated in 
any part of the system.* It can be 
shown that the number of w ays a 
retention time “R” can be obtained by 
summing “K” components pj, P., 
---- px Where all R, pj, are integers is 
given by: 


KR = 1 (R - 1)! 
K -— 1 


(K —- (R — K)! 
Thus the probability of a particle 
remaining in a “K” agitator system a 


total time “R” is given by: 


R-1 
Pr(R) = (j rK¥ (1 — r)R-K 
K -1 


Alternatively of an original “n” 
particles entering the first agitator at 
time zero, the number leaving during 
the Rth time period will be: 


R-1 R 1 
y = ). are — r)R-K.,......(2) 
R K — 1 


This equation also gives the number 
of particles leaving the Kth agitator 
per unit time which have resided i 
the system a total time “R” 


It should be noted that this result 
is valid for any agitator in the sys- 


_°We are assuming here that time intervals 
hdve been chosen sufficiently small that the 
retention time of a particle in any one agitator 
will be an integral number of time units (i.e. 
no fractional times). Thus a particle cannot be 
retained in any agitator less than one time unit, 
and hence in a “K” agitator system, the mini- 
mum retention time must be “K” time units 
which can only be obtained in one way, i.e. one 
time unit in each agitator. It there fore follows 
that for the particle retention time “R” to be 
obtainable in more than one way, “‘R” must be 


greater than K. (R>K). 
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Figure 


tem, not only for the last. By a simple 
transformation it can be shown that 
this is a Poisson distribution. 


Figure 2 illustrates this function for 
n= 100,000, r= 0.01 and K = fF, 2, 3, 
4, and 5. 


This distribution function — re- 
sembles results obtained by others 
considering the purging of a series of 
stirred tank reactors. 


2. Growth Function for Crystals 
an Agitator System. 
The following development of a 
growth function is similar to that 
used by W. C. Saeman ‘!), 


Assume that all agitators in a sys- 
tem are held at a steady degree of 
supersaturation “S” and ‘further that 
“S” is held at such a level that al- 
though crystals in the suspension will 
grow, no nucleation will take place. 
It is assumed that nuclei are provided 
by seeding in the first agitator. 

Let the volume of a single crystal 
be v=k,L? where “L” is a linear 
dimension of the crystal and “k,” is a 
shape coefficient. 

Assume that the linear rate of 
growth of a crystal is proportional to 
the supersaturation of the solution in 
contact with it, Le.: 

dL 


k2S 
dt 


where k, is a coefficient relating rate 
of growth to supersaturation, In- 
tegrating and assuming that the origi- 
crystal 


nal size of the nuclei is 


negligible, w 


e have 


L = k St (3) 
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Therefore the volume of a crystal 
after growing a time “t’’ is: 

ve = ki (k2St)3 (4) 
Further if we let the density of the 
crystal be “p” the weight of a crystal 
that has been growing for time “e” is: 


Wt = pk, (k2St)? (5 
Development of the Size Distribu- 
tion Function for Crystals Grown 
in a Series of Agitators. 

(a) General 
The results of sections 1 and 2 will 

now be combined retaining all the 

previous assumptions. In particular it 

is assumed that the supersaturation “S” 

is maintained in each agitator at the 

same level such that crystal growth 
proceeds but no nucleation occurs. 

It is appreciated that in practice a 

constant inflow and outflow of “u” 

per unit time for each agitator would 

be difficult to achieve exactly while at 
the same time maintaining “S” con- 
stant. Such a condition might be ap- 
proached either by cooling the solu- 
tion or by evaporating quantities of 
solvent equal in volume to the solute 
input to each agitator. We will as- 
sume that this condition is achievable. 


(b) Total Weight and Size Distri- 
bution of Crystals Discharged 
from the Kth Agitator 


Ir has been shown that the number 
of crystals leaving the Kth agitator 
which have been retained in the sys- 
tem a time “R” is given by equation 
(2). The weight of a crystal which 
retained in the system a 
is shown by equation (5) 


has been 
time “R” 
to be: 

Wr = pk, (k2SR)5 
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TABLE 1 


K 


VALUES OF Z 


1 2 
a ib Z 
a 1 
0.5 99 .82 99.98 
1 98.1 99 63 
2 85.7 94.7 
3 64.7 81.5 
4 | 43.4 62.9 
5 26.5 44.1 
6 15.1 28.5 
7 8.2 17.3 
8 4.2 10.0 
9 2.1 5. 
10 1.00 2.3 
11 OB | 
12 0.76 
13 
14 
15 
16 
Combining these two equations, 


the weight of crystals leaving the Kth 
agitator w hich have been retained in 
the system a time “R” is given by: 


R-1!1 
pk; (k2SR)*. 
K —! 


x nrk (1 — r)R-*K 


Wr = 


The total weight of crystals leav- 
ing the system per unit time which 
have been retained in the system a 
time “R,” or less is: 


W61_= pk: (k2S)* nr¥ 7 es 


R -1 
. R?(1 — r)®-* (6) 
K — 1 


If the time units are made suff- 
ciently small this equation may be ex- 
pressed as an integral. Making the 
time units very small also makes “R” 
large for all crystal sizes of practical 
interest and at the same time makes 
“r” very small. The retention time of 
any crystal may be expressed as a 
multiple of the mean retention time 
per agitator, i.e. 


R =a/r where0 Sa <+ om 


When ‘R” is large: 


R-1\). RE — (a/r)K-! 
cat? CK-m ee — 


and 


(1 — r)®-* = (1 —a/R)®-* = e- 


Equation (6) can therefore be writ- 
ten: 
WE: = W2! = pki (k2S)? nrX. 

(a/r)=—! da 


i (K — pr 62/0? Ce) 
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a 


3 4 5 

Z Z Z 
1 a a 
99 94 99 99 99.99 
98.3 99.55 99.89 
91.6 96.65 98.8 
78.5 88.93 94.9 
61.6 76.22 86.7 
44.6 60.6 74.4 
30.1 45.0 59.9 
19.1 31.3 45.3 
11.6 20.7 32.4 
6.7 13.0 22.0 
3.8 7.9 14.3 
2.0 4.6 9.0 
1.07 2.6 5.4 
1.42 $2 
0.76 1.80 
1.00 


. da 


pk, (keS)* n 
(K — 1)! r 


where H = 


Integrating it is found that the 
weight of crystals leaving the Kth 
agitator per unit time which have 
been retained in the system “a” 
mean retention times or less is given 
by: 


Ww: 
| Ke 2)t =e Z, Pi aero] 
(7) 
ital p* " 
(K + 2) (K + 1) (K +2 —i+1) 


is the number of permutations of 
(K +2) things taken “i” at a time. 

For the special case of K = 1 this 
same result has been obtained by 
Saeman (1), 

The total weight of crystals leaving 
the Kth agitator per unit time is ob- 
tained when we allow “a” to approach 


infinity in equation (7) 


ce 


W, = H(K +2)! 


°o 


The fraction F, of the system dis- 
charge per unit time that has resided 
in the system for “a” mean retention 
times or less is of course given by: 


w: /w® 


and this fraction times 100 gives the 
complement of the usual cumulative 


Fa 


weight % figure used in describing a 
size analysis. Therefore letting Z* 
represent the percentage by weight of 
effluent crystals which have resided in 

‘0 


the K agitator system at least “a 
mean retention times: 


ZK = 100(1 — Fy = 


100.e-* St? kt? 
4 (K+2—i) 
(K + 2)! a Pia om 


This expression gives the cumula- 
tive weight % of crystals in the dis- 
charge of the Kth agitator which are 
of a size greater than “L” where: 


‘ = koSa r 


The important fact shown by these 
equations is that the cumulative 
percent w eight of crystals which 
have grown in the system for “a” 
mean retention times or longer de- 
pends only on the number of agita- 
tors. The size of the crystals will of 
course depend on the length of the 
mean retention time and the crystal 
growth rate (k,S). Cumulative dis- 
tribution curves for K = 1---5 are 
shown in Figure 3 and values of Z* 
for K = 1---5 and a = 1---16 are given 
in Table 1. 


It can easily be shown that the 
points of inflection of the curves 
(which correspond to the maxima of 
the non-cumulative distribution) oc- 
cur when a= K+2. This implies that 
the greater the number of agitators 
in the system the closer the maximum 
approaches to the overall mean reten- 
tion time of the system. 


It is interesting to note that in a 
single agitator 50% by weight of the 
crystals have been retained at least 
3.7 times the overall mean retention 
time. In a series of 5 agitators, 50% 
by weight of the crystals have been 
retained at least 1.54 times the overall 
mean retention time of the system. 
This implies that for equiv alent sys- 
tem size the single agitator will grow 
considerably larger crystals than a 
series of agitators. 


(c) Calculation of Other Functions 
of Interest 


Using a development similar to that 
set forth above it is easy to derive the 
following expressions: 

(1) Number of nuclei per unit 
weight of crystals leaving the Kth 
agitator: 

_ r?.(K — 1)! 
pki (ksS)® (IK +2)! 


fase 
N/ Wtot 


(ii) Gain in crystal weight in unit 
time in the Kth agitator: 


3pk; (k2S)% nk (K + 1 


r3 


AW = 
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From this it can readily be shown 
that for every unit gain in weight in 
the first agitator of a series, the gain 
in the Kth agitator is K(K+1). 


2 


(ii) Surface area per unit weight 
of crystals leaving the Kth agitator: 


= 


A/W = 
(K + 2). pkik:S 


where k, is a shape constant such that 
for one crystal A = k,L?* 


4. Comparison of Observed and Cal- 
culated Size Distributions. 

At the Trail operations of Cominco 
ammonium sulphate for fertilizer use 
is produced by evaporation of solu- 
tion in forced circulation vacuum 
crystallizers. If we regard such an 
evaporator as a single agitator in 
which perfect mixing takes place we 
have from equation (8) that the 
cumulative percent weight of crystals 
in the discharge which have been 
resident a time a/r or longer is given 
by: 

100e-* 


Zi = (a8 


31 + 3a? + 6a + 6) (9) 


In the course of a series of tests a 
large number of screen analyses of salt 
sampled at the cry stallizer discharge 
were obtained. After eliminating 
those taken on days earlier than the 
3rd after start-up, to ensure that 
equilibrium had been reached, 61 sets 
remained which gave the average 
screen anlaysis shown in Table 2. 

It is not a difficult matter to calcu- 
late the theoretical distribution ac- 
cording to equation (9) which mini- 
mizes the sum of squared deviations 
of the experimental data. 


A sufficiently good approximation 
can usually be more readily obtained 
by drawing a smooth curve through 
the data, choosing a point on this 
corresponding to the percentage point 
for some value of “a” which is an 
integer, and reading the remaining 
points off from Table 1, using these 
to draw the theoretical curve. 


The actual and best fitting theo- 
retical distributions are presented i 
Table 2 and are _ illustrated in 
Figure 4. 

At one time ammonium sulphate 
was produced at Cominco’s Trail 
plants by direct ammoniation of sul- 
phuric acid in saturators, which were 
merely single stirred agitators. The 
comparison between actual and theo- 
retical size distributions for these 
units is also shown in Table 2 and 
Figure 4. 


TABLE 2 
COMPARISON OF ACTUAL AND THEORETICAL SCREEN ANALYSES 


VACUUM CRYSTALLIZER 


Screen Size (Tyler) +14 


Cum. % W t. 
Actual 
Theoretical 


60.4 
60.9 


82 
80 


+20 +28 +35 +48 


93.6 97.4 98.7 
92.1 bY 99.1 


sun 


SATURATOR 


Screen Size (Tyler) +20 


Cum. % Wt. 
Actual 14.6 
Theoretical 16.1 


99.99 | 


999 


99.8 
99.5 
99 


98 






95 


90 


80 


7O J 
60 
50 
40 


30 | 


20 4 


CUMULATIVE PERCENT WEIGHT 


CALCULATED POINTS 


CALCULATED POINTS 


FOR SATURATOR 
o2 


0.1 


wel 


@ 
nN 





+ 100 MESH 
35 


- 20 


4 7 — 


+ es 
r 48 


a GED. 


CRYSTAL 


VACUUM CRYSTALLIZER 


41.9 68.9 
41.0 66.0 


— OBSERVED DISTRIBUTION FROM 





OBSERVED DISTRIBUTION 


FOR SATURATOR 


FOR VACUUM CRYSTALLIZER 


10 MESH 


< 


yy 
° eo) 0.02 0.03 0.04 


00s 0.06 0.07 


Size — InNcHES 


Figure 4 


These two sets of comparisons 
show that the theoretical develop- 
ment gives excellent agreement with 
the results obtained in practice on 
two units of widely different con- 
struction, size and operating condi- 
tions, operated as single vessels. No 
experimental data are available for a 
cascade of agitators, but since the 
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theoretical development fits the single 
agitator data so well, there is every 
reason to believe it will be valid for 
agitators operated in series. 


5. Some Comments on Nucleation 


and Seeding. 


In the derivation of the preceding 
equations it was assumed that there 
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would be a constant rate of seeding 
into the first agitator of a series of 
“n” nuclei per unit time of negligible 
size. In pri actice such a procedure is 
seldom, if ever, followed. It was cer- 
tainly not applied in the operating 
units which were 


used for compari- 


son with the theoretical size distri- 


bution function. Nevertheless “n” can 
have an indirect effect in that it can 
influence the level of “S” and hence 


the growth rate. 


It is generally agreed that below 
some level of supersaturation, “S,”, 
little if any Above 


this level the rate of nuclei formation 


nucleation occurs. 


increases very quickly with increase 
in supersaturation. An 
the type: 
n= b(S 


equation of 


S:)* where S2S, and x>1 


110 


5 


might perhaps be used to describe the 
rate of nuclei formation. Some rather 
fragmentary eee work at 
Cominco indicates that “x” has a value 
greater than 4 in the case of am- 
monium sulphate, and that in the pre- 
of normal turbulence such as 
would be encountered in any 
industrial crystallizing unit it is 
difficult to ratio 


sence 


very exceed the 
s-S is 

——— = 5/4 where S, is a saturated 
of. 


solution. 

The net effect of this behaviour is 
that in a crystallizing unit in the 
absence of deliberate seeding, the level 
of “S” tends to be always just a little 
above “S,” and is very stable for any 
given set of operating conditions. For 


if “S” falls below this level there will 
be a marked decrease in the rate of 
nucleation and the surface area in the 
unit will decrease until “S” rises to 
the stable level, and similarly if “S” 
increases above this level the rate a 
nuclei formation will increase very 
rapidly which will increase the crystal 
surface and eventually result in a 
decrease of “S” to the stable level. 

Therefore in a single agitator sys- 
tem “S” will probably change very 
little over a wide range of operating 
rates. In a cascade of agitators, the 
first one of the series will tend to 
produce nuclei at a stable rate, thus 
giving the equivalent result to seed- 
ing, but care will have to be taken 
that the supersaturation in the suc- 
ceeding agitators remains below the 
level where nuclei formation occurs 
or else the preceding mathematical 
model will not be applicable. 


Effect of Production Rate on Crys- 

tal Size Distribution. 

The effect of a change in produc- 
tion rate of a crystallizing unit on the 
crystal size distribution can be de- 
termined readily with the mathe- 
matical model which has been de- 
veloped. The major direct effect of 
a change in operating rate is to alter 
the mean retention time in the system 
by a corresponding factor. This is of 
course assuming that percentage of 
solids in the magma and the total 
W eight of material in the system re- 
main practically unchanged, which is 
generally the case in an_ industrial 
unit. If our argument in the preced- 
ing section is valid, then the growth 
rate (k,S) may be expected to change 
only very slightly when the produc- 
tion rate is changed. Therefore, since 
the mean retention time varies in- 
versely with the production rate, and 
the crystal growth rate remains rela- 
tively constant, the crystal size at- 
tained in any given multiple “a” OF 
the mean retention time will be in- 
versely proportional to the unit pro- 
duction rate. Hence if, at production 
rate “X,” we have available a graph 
of crystal size against cumulative 
weight percent, to find the size dis- 
tribution at production rate “X,,” the 
crystal size corresponding to any 
cumulative w eight percent need only 


be multiplied by ‘ "ay / ain « 
During the period when am- 


monium sulphate was produced at 
Cominco’s Trail operations in both 
vacuum crystallizers and saturators, 
the crystallizers produced a 
siderably larger product than the 
saturators. The saturators however 
were operated at more than twice the 
production rate per unit volume 


con- 
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compared with the crystallizers. It is 
interesting to see how the product 
size from the two units would have 
compared if both had been operated 
at the same rate per unit volume. The 
saturators were operated at a rate of 
78 tons per day, and to be equivalent 
to the Crystallizer operating rate, 
this would have to be reduced to 33 
tons per day. Figure 5 shows the 
crystal size distribution for the 
vacuum crystallizers and for the 
saturators as operated at 78 tons per 
day. Using the method described 
above the expected distribution at 33 
tons per day is shown on the same 
figure as a broken line. It is apparent 
that at similar production rates per 
unit volume the two types of unit 
would be expected to produce 
crystals of comparable sizing. 


Prediction of the Performance of a 

§-Agitator Crystallizing System. 

This study was initiated in an at- 
tempt to predict whether an existing 
system of agitators should be operated 
in series or in parallel as saturators for 
ammonium sulphate production. The 
result of our study in this connection 
is presented in this section. 


The system consists of a cascade of 
5 similar agitators each with a volume 
of 1600 cubic feet. It will be assumed 
that we wish to produce a product 
which contains 60% by weight of 
+ 14 mesh material. It is required to 
know the tonnage of this material 
which can be produced operating the 
system in series and in parallel and 
also. the complete size distribution 
for each case. 


(a) Estimation of Constants 


p = density of ammonium sulphate 


110 Ibs per cu ft 
k, = shape constant relating linear di- 
mension and volume 0.866 
keS = rate of growth of ammonium sul- 
phate per unit time. 

In the analysis, k, and S only occur 
in the form of their product and the 
value of this can be estimated from 
available operating data for saturators 
and crystallizers. 

From the volume of the saturators 
previously operated, the production 
rate and the magma density, it was 
estimated that the 


time in the saturators at a 


retention 
rate of 78 
tons per day was 293 minutes. From 


mean 


lable 1 it is found that in a single 
64.7% by weight of the 
erystals will have resided in the sys- 


agitator 


tem at least three times the mean re- 
tention time. From Figure 4 it is 
found that 64.7% by weight of the 
crystals have a linear dimension of at 
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least 0.0175 inches. 
tion (3) 


Thus from equa- 


0.0175 
3 xX 293 «* 12 


kS = 


1.66 X 10-° feet per minute 


Similar calculations for the vacuum 
crystallizer give a value of: 


kS 2.00 10-° feet per minute 


Since the two types of unit are very 
different and also the determination of 
magma density is not very accurate, 
a difference of this amount is not sur- 
prising. 


(b) Operating rate to produce 
60% +14 mesh with series 
arrangement 

From Table 1 it is seen that for a 

cascade of five agitators the 60% 

point in the cumulative weight distri- 

bution function occurs at about 7 

mean retention times. Therefore a 14 


mesh particle (0.046 inches) must 
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grow in 7 mean retention times. The 
estimated value of k,S for a saturator 
is 1.66 x 10~® feet per minute. Then 
by equation (3) the mean retention 
time must be: 


0.046 
4 M12 XM LOO XM IO 
330 minutes 
Usually in a crystallizing system 
the slurry containing the crystals 
must be pumped from the last agita- 
tor to centrifuges or filters, so the 
maximum concentration of crystals in 
the last agitator is set by the limita- 
tions of the pump. In Cominco prac- 
concentration of 42 Ibs. of 


ammonium sulphate crystals per cubic 


tice a 


foot of slurry has been found to be a 
reasonable operating maximum. Using 
this figure, the total weight of crystals 
suspended in the last agitator is: 


1600 « 42 67,200 Ibs 


Ii! 


TABLE 3 


COMPARISON OF SIZE DISTRIBUTIONS OBTAINED IN SERIES AND PARALLEL OPERATION 
OF A 5-AGITATOR SYSTEM 


Screen Size 


(Tyler) +10 +14 +20 
Com. % Wt. 
Series 147 
T/D 23.7 60.0 87.5 
Parallel 335 
T/D 34.0 60.0 80.0 


Mean retention time = 
Weight of suspension 


Production rate 


67,200 


Therefore duction rate = = 
Therefore productio 330 


204 Ib/min = 147 tons/day 


Thus with 5 agitators, each of 1600 
cu. ft. capacity, operated i in a series, a 
production rate of 147 tons per day 
of ammonium sulphate could be ex- 
pected having a screen analysis of 
60% + 14 mesh. 


From equation (8) the screen size 
distribution curve is readily cal- 
culated, and this is illustrated in 
Figure (6). 

The input of raw materials, am- 
monia and sulphuric acid, to each 
agitator can be easily calculated from 
the equation presented in section 
3-c-ii. 


(ce) Operating rate to produce 
60% + 14 mesh with parallel 
arrangement 

Proceeding as in the previous sec- 
tion, it is found that the 60% point of 
the cumulative weight distribution 
curve for a single agitator system 
occurs at 3.2 mean retention times. 

Therefore the mean retention time in 

an agitator must be: 


0.046 
re 12 x 3.2 X 1.66 X 10-* | 


- 


722 minutes 


The production rate for a single agi- 
tator is: 


67,200 
a= = 93 lb/min 
722 


= 67 tons per day. 


Therefore for 5 similar 


operated in parallel: 


agitators 


112 


+28 +35 +48 + 100 
97.4 99.55 99.93 99 99 
92 3 97.3 99 20 99.91 
Overall production rate = 335 tons 


per day. 

As before the screen size distribu- 
tion may be calculated by equation 
(8) and it also is shown in Figure 6. 
The comparison of the screen size 
distribution of the two methods of 
operation is also given in Table 3. 

It will be noted that when the pro- 
duct contains 60% + 14 mesh material 
the series system produces slightly 
less fine material, but the penalty in 
reduced production rate is very 
heavy. In fact it can be shown that 
at similar operating rates for the two 
systems, both will produce a similar 
amount of fines but the parallel opera- 
tion will produce a much higher per- 
centage of the large screen sizes. 

It must be concluded therefore that 
series operation of a number of crys- 
tallizing units is not practical for pro- 
duction of fertilizer grade ammonium 
sulphate. 


Conclusions 

In a crystallizing system consisting 
of one or more well mixed vessels, 
the crystal size distribution is de- 
scribed by a simple function whose 
only variable is the mean retention 
time. 

The mean residence time of crystals 
in a growing system is much greater 
than the mean residence time of the 
overall material in the system. For 
example 50% by weight of the 
crystals produced in a single agitator 
have stayed in the system at least 3.7 
mean retention times. 

For a given percentage of large 
sized crystals, a series of vessels will 
produce a little less fine material than 
a single vessel of the same total 
volume, but the output from the 


series arrangement will be very much 
less than from the single vessel. 


Suggestions for Further Study 


The authors feel that the approach 
used in this study could be extended 
to more complicated crystallizing 
problems. A few of the problems 
which come to mind are: 


(a) The effect of a magma recycle in 
a multi-agitator system. 

(b) The effect of settling and pre- 
ferential removal of large crys- 
tals. i ; 

(c) The effect of zones of varying 
supersaturation within a system 
which would lead to non-uniform 
growth or even to disappearance 
of crystals: this might perhaps be 
handled by a stochastic function. 


(d) The effect of transient distur- 
bances to the system. 


Some of these problems would un- 
doubtedly lead to much more com- 
plex models, but this is not an in- 
superable difficulty if access could be 
had to a large digital computer. 


Nomenclature 

a = number of mean retention times. 

e = base of natural logarithms. 

k = constant. 

K = number of agitators in a cascade. 

L = linear crystal size. 

n = number of particles entering and 
leaving a cascade per unit time 
under equilibrium conditions. 

p = density of crystal. 

R = retention time of a particle in a 
cascade, 


1/r = mean retention time in one agitator 
of a cascade. 

S = supersaturation. 

volume of material entering and 

leaving a single agitator per unit 

time. 

volume of a single agitator. 


vy = volume of a single crystal. 

w = weight of a single crystal. 

W = weight of crystals leaving a cascade 
per unit time. 
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Circulation and End Effects 


ina 


Liquid Extraction Spray Column 


D. CAVERS?® and J. E. EWANCHYNA®* 


The operation of a spray type liquid-liquid extraction column was 
studied. Acetic acid was transferred between a continuous water phase 
and a dispersed methyl isobutyl ketone phase. 


Both the continuous phase and the dispersed phase were sampled at 
elevations between the spray nozzle at the bottom of the column and the 
interface at the top. Thus concentration profiles for each phase were 
obtained. 


Discontinuities in the concentration profiles were noted at the 
interface. For transfer of acetic acid from the aqueous phase to the 
ketone phase the discontinuity in the water phase concentration profile 
could be broken into two parts: one representing the effect of drop 
agitation at the interface, and the other the effect of back-mixing in the 
aqueous phase. The discontinuity in the ketone phase concentration 
profile was due only to the agitation effect. For transfer of acetic acid 
in the reverse direction the discontinuity in the ketone phase profile 
was absent, and that in the water phase was due only to back-mixing in 
the continuous phase. 


The concentration profiles were used to calculate overall capacity 
coefficients and H.T.U.’s. 
Introduction a technique for sampling the continu- 
ous phase of a spray extraction column. 
From the data obtained. they calcu- 
lated “internal” mass transfer coeffi- 
cients for short lengths of column. 
Plots of concentration of the continu- 
ous phase versus distance from the 
spray nozzle made by these workers 
showed little or no inlet end effect 
for the entering d’ spersed phase al- 
though a large effect of this sort might 
have been expected from results such 
as those of Sherwood, Evans and 
Longcor “), and West, Robinson, 
Morgenthaler, Beck, and McGregor‘®), 
For the most part the points in each 
of these concentration plots formed a 
smooth curve including the inlet con- 
centration of the drops as one point, 
no abrupt changes in concentration or 
only minor ones being observed as the 
drops were formed. However, a con- 


IQUID-LIQUID extraction is a chemical 
engineering unit operation which 
has been used in industry with in- 
creasing frequency in recent years. 
The increased importance of liquid- 
liquid extraction has been reflected in 
a considerably increased interest in 
research in the various aspects of ma- 
terial transfer in liquid-liquid systems. 
A fairly large portion of this research 
has dealt with spray liquid extraction 
columns because of their inherent 
simplicity. 

Mass transfer coefficients were de- 
termined in most cases by using the 
inlet and outlet concentrations and 
calculating the logarithmic mean driv- 
ing force. However, Geankoplis and 
Hixson “); Geankoplis, Wells, and 
Hawk (); Kreager and Geankoplis “); 
and Vogt and Geankoplis ‘) perfected 
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siderable end effect was found at the 
continuous phase inlet to the column. 
An unexpectedly large portion of the 
total extraction took place when the 
drops of dispersed phase coalesced at 
the interface, and, as a result, the 
continuous phase concentration profile 
showed a discontinuity at this point. 
Geankoplis and Hixson “) attributed 
this end effect to turbulence caused 
by the coalescence of the drops, but 
suggested that changes in interfacial 
tension might also be important. 


Keith and Hixson ‘?) explained the 
absence of an end effect at the dis- 
persed phase inlet in work such as 
that of Geankoplis and Hixson!) as 
being due to a rapid rate of drop 
formation and a consequent short resi- 
dence time for each drop at the 
nozzle. 

Newman (8) suggested that the end 
effect found by Geankoplis et al. (1.2) 
at the continuous phase inlet was the 
result of vertical mixing of the con- 
tinuous phase due to the movement 
of the drops. Newman showed that 
some of the data of Geankoplis et al. 
was midway between the results to be 
expected for truly countercurrent, un- 
mixed flow of the continuous phase, 
and the completely uniform concen- 
tration in that phase which would be 
produced by very efficient stirring. 
Although considerable research has 
been done on circulation within drops, 
mainly by Garner and his co-workers, 
e.g. as described in References (9),(10), 
(11),(12), and (13), little work has been 
done on circulation in the continuous 
phase. 

However, this problem was investi- 
gated in a 6 in. inside diameter liquid 
extraction column by Gier and 
Hougen (4), These workers sampled 
both the continuous and dispersed 

hases at various locations. Their re- 
sults showed considerable back-mix- 
ing of the continuous phase when the 
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column was operated as a spray tower. 
When it was operated as a packed 
tower the mixing was reduced. Gier 
and Hougen pointed out that in deter- 
mining the overall height of a transfer 
unit, ‘graphical integration is neces- 
sary, this integration being dependent 
ona knowledge of the concentration 
of each phase throughout the column. 


Heertjes, Holve, and Talsma (!5) 
measured the concentration of the 
continuous phase of a spray column 
by draining the column at the end of 
runs and sampling the effluents. Dia- 
meter to height ratios in the range 
0.057 to 7.3 were investigated. The 
concentration of the continuous phase 


was nearly constant over the column 
height, indicating good longitudinal 
mixing. These results do not agree 
with those of Geankoplis et al (12:3), 

The present investigation was begun 
in an attempt to shed more light on 
circulation in the continuous phase of 
spray columns and on end effects. 
Less circulation might be expected in 
columns of lower diameter to height 
ratio (14), The column described in the 
present paper had an inside diameter 
of 1.5 in. The diameter to height ratio 
was approximately one third of that 
used by Gier and Hougen. A method 
for sampling each of the phases at 
various elevations in the column was 





Figure 1—Schematic flow diagram. 
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developed from the methods used by 
Geankoplis and Hixson (!) and Gier 
and Hougen ‘'#), 


The system investigated was water- 
acetic acid- methyl isobutyl ketone. 
The direction of transfer of acetic 
acid was usually from the water to the 
ketone. Experimental runs in which 
this direction of transfer applied have 
been called normal runs in this paper. 
Those in which the acetic acid was 
transferred from the ketone to the 
water have been called reverse runs. 


Experimental 


Materials. The methyl isobutyl 
ketone used was technical grade ob- 
tained from Carbide Chemicals Com- 
pany, Toronto. Glacial acetic acid, 
C.P., was obtained from Canadian 
Industries Limited, Toronto. The 
water used was laboratory distilled. 


Apparatus. Figure | is a schema- 
tic flow diagram of the spray-type 
liquid extraction column used in this 
research, together with accessories for 
maintaining, controlling and measur- 
ing the flows of the two phases, and 
for sampling them. The methyl 
isobutyl ketone drops were sprayed 
into the continuous aqueous phase at 
the bottom of the column. The con- 
tinuous phase entered the column at 
the top and flowed downward, 
countercurrent to the rising ketone 
drops. Except as noted in the follow- 
ing discussion the flows were carried 
in’ Saran tubing, in stainless steel 
Schedule 40 pipe, or in Pyrex brand 
“double tough” pipe. All stainless 
steel used was Type 304 and all gas- 
kets were asbestos. 


In Figure 1, F and S were feed 
tanks situated on the floor of the 
laboratory and used to store the 
aqueous phase and the ketone phase 
respectively. A Sigmamotor tubing 
pump P was used to elevate the phases 
to small constant head tanks A and B 
made from stainless steel. The over- 
flows from these tanks were returned 
to the respective feed tanks. The 
water phase flowed by gravity from 
constant head tank A through cali- 
brated rotameter D and “entered the 
top of the column at N through two 
4 in. stainless — steel pipes. The 
ketone phase flowed by gravity from 
constant head tank B through cali- 
brated rotameter G, and into the 
column through spray nozzle T. 
Stainless steel needle valves, K and M, 
were used to control the rates of flow 
of water and of ketone respectively. 
After entering the column at N the 
water phase flowed downward 
through the column R, out at the 
bottom, and into storage tank U. 
After being formed at spray nozzle 
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T the ketone drops rose through 
column R_ countercurrent to the 
descending water phase. The drops 
coalesced at interface Q to a con- 
tinuous methyl isobutyl ketone phase 
which left the column at the top and 
flowed into a storage vessel E. Ves- 
sels F, S, U, and E were all twelve 
U.S. gallon capacity Pyrex carboys. 
Vent to atmosphere W on the ketone 
discharge line was used to prevent 
siphoning and resultant irregular 
flows of ketone out of the top of the 
column. The position of interface Q 
was controlled by varying the height 
of vent X on the water phase dis- 
charge line. 

Figure 2 shows the column proper 
in more detail than does Figure 1. 
The design is similar to that of 
Blanding and Elgin “'%). A standard 
Pyrex reducer was used up to the end 
of Run 34 to provide space in the 
column for the spray nozzle. For 
Runs 35 to 38 the standard reducer 
was replaced by the conical connector 
shown in Figure 3. For Runs 39 to 44 
inclusive the reducer and the two 
lengths of 14 in. pipe above it 
(Figure 2) were replaced by the one- 
piece column section also shown in 
Figure 3. The standard reducer had 
the disadvantage that the change in 
diameter from 13 in. to 3 in. took 
place rather suddenly, whereas the 
change was much more gradual in the 
devices shown in Figure 3. The top 
of the nozzle was so situated with 
respect to the reducer in use that the 
free annular area around the nozzle 
was always at least equal to the cross- 
sectional area of the column _pro- 
per (3), The concentration profiles ob- 
tained in a run when the standard re- 
ducer was used were compared with 
concentration profiles for the same 
run repeated after the installation of 
the connector shown in Figure 3. No 
appreciable change was observed in 
the profiles. The upper end of the 
1} in. Pyrex pipe forming the top of 
the column was plain: the usual in- 
crease in outside diameter necessary 
for attaching a Corning flange was 
absent. 


Two nozzles were used in the pre- 
sent study for dispersing the ketone 
drops. These are shown in Figure 4 
and Figure 5 respectively. The nozzle 
design was based on that of Kreager 
and Geankoplis (8), Nozzle I being 
identical in most respects with theirs. 
The nozzle tips were made from § in. 
O.D. by 0.020 in. wall thickness stain- 
less steel tubing. The inside diameter 
was increased to 0.100 in. by drilling. 
For the present Nozzle Il, the tips 
were chamfered at 45° to the sharp 
edge recommended by Hayworth and 
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Figure 2—Detailed drawing of column. 


Treybal "7, The flow of liquid was 
divided approximately evenly among 
the several tips. In the case of Nozzle 
Il, this distribution was checked by 
counting drops. Nozzle | was used to 
the end of Run 18 and Nozzle Il for 
Runs 19 to 44 inclusive. 

With the object of obtaining drops 
of uniform size, Hayworth and 
Treybal '7) recommended that the 
linear velocity of the fluid in the 
nozzle tips should be held constant at 
a value below approximately 0.3 
ft./sec. In the present work this 
velocity was maintained in the range 
0.266 to 0.296 ft./sec. A velocity in 
this range was used also by Kreager 
and Geankoplis “). Continuous 
streams of the dispersed phase liquid 
approximately 0.75 in. long existed in 
the present study between the nozzle 
tips and the respective points of drop 
detachment. Thus the ketone velocity 
was above the “jetting point” ‘7. 
Different ketone flow rates were ob- 
tained by making use of different 
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numbers of nozzle tips. Those not in 
use were blocked off by means of 
small corks, or later, by means of 
Teflon caps. Figure 6 shows these caps 
and also the patterns used for block- 
ing nozzle tips at various ketone flow 
rates through Nozzle II. A run made 
using Nozzle I was repeated using 
Nozzle II. Little or no change in the 
concentration profiles resulted. 

When the column was in operation 
at the steady state, samples of each 
of the phases were withdrawn from 
various elevations above the nozzle 
by sample probes shown at b and a in 
Figure 1. A partial vacuum was used 
to withdraw the samples. Two dif- 
ferent sets of sampling tubes were 
used. 

The first set was made from § in. 
O.D. and 0.08 in. 1D. Pyrex tubing, 
and the second set from in. O.D. by 
0.085 in. 1.D. stainless steel tubing as 
shown in Figures 7 and 10. The Py rex 
tubes were of similar design. ‘To- 
gether the Pyrex hook and inverted 
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funnel ocupied 4.2% of the internal 
cross-sectional area of the column. 
The corresponding figure for the 
stainless steel sampling tubes was 
6.2% 


The sample tubes passed through a 
Saran sleeve in the top end plate of 
the column. (See Figure 2.) They 
were attached to a w ooden block 
shown at c in Figure 1. The elevation 
of the sampling tubes in the column 
was changed by raising or lowering 
the block in guide d. A bolt through 
the block and the guide held the ends 
of the tubes at certain definite dis- 
tances above the nozzle tips. 

The sampling tubes were attached 
to 250 ml. erlenmeyer flasks at h and 
k in Figure 1 by means of & in. LD. 
by 3/64 in. wall thickness rubber 


PYREX BRAND 
"DOUBLE TOUGH" 
10° REDUCER 


tubing. A vacuum of 2.5 to 3.5 in. of 
mercury produced by a water aspira- 
tor served to draw the two samples 
into the erlenmeyer flasks. The 
sampling rates were controlled by ad- 
justing screw clamps e and f in 
Figure 1. 

The temperature of each of the 
inlet and the outlet phases was 
measured in the early runs by thermo- 
meters suspended in the liquid i in the 

carboys at F, S, U, and E in Figure 

1. In the later runs, these tempera- 
tures were measured by means of 
Yel-O-Bak midget industrial ther- 
mometers with 18-8 stainless steel 
wells in direct contact with the 
fluids. These thermometers are shown 
at J,, J.. J,, and J, in Figure 1. 

The use of the Sigmamotor tubing 


|. 4.07"__»| 


PYREX BRAND "DOUBLE 
TOUGH" 


COLUMN WITH 
10° REDUCER 





Figure 3—Tapered connector and one-piece column. 
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pump referred to previously, and lo- 
cated at P in Figure 1, requires that 
each of the phases pass through a 
length of flexible tubing. Sufficiently 
flexible tubing not attacked by methyl 
isobutyl ketone was unav ailable. Red 
natural rubber tubing, and, later, 
neoprene tubing were used in the 
pump. Both contaminated the ketone 
to some extent, the natural rubber 
probably more so than the neoprene. 
The extent of contamination in- 
creased from run to run as the ketone 
was reused. However, the impurities 
did not appear to affect the extraction 
data. 


Procedure. Johnson and Bliss (18) 
showed the existence of a small effect 
of solute concentration on the overall 
mass transfer coefficient. In the pre- 
sent work the solute concentration of 
the feed solutions was kept roughly 
the same for each of the runs so that 
any effects of concentration should 
be negligible. 


To insure that the water phase was 
saturated with ketone a thin layer of 
the latter was maintained in the water 
phase feed carboy. The layers were 
stirred together for at least two hours 
and were allowed to settle for at least 
one hour before a run was begun. 
Similar precautions were taken with 
the ketone phase. 


At the start of a run, ketone was 
allowed to flow slowly through the 
nozzle while the column was filled 
with aqueous feed to the point where 
the interface was to be maintained. 
This procedure prevented the water 
phase from entering the nozzle. Both 
flows were then set at the desired 
rates by means of needle valves at K 
and M in Figure 1. The height of the 
interface was adjusted by either 
raising or lowering the vent at X in 
Figure 1. 

Steady state conditions were 
reached in the column after sufficient 
of the feeds had been added to change 
the contents of the tower two to 
three times. Internal samples were 
taken through the sample tubes des- 
cribed previously. The rate of samp- 
ling was 10 ml./min. and the duration 
approximately seven minutes. Experi- 
ments showed that four minutes 
purging was sufficient to change the 
contents of the lines so that the 
material collected in the first four 
minutes was discarded, and that col- 
lected in the last three minutes re- 
tained as the sample. Effluent phase 
samples were taken from time to time 
while internal sampling was in pro- 
gress. 


The velocity of the rising ketone 
drops was determined several times 
during each run by measuring the 
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time required for a ketone drop to 
pass between two marks on the side 
of the column. 


The amount of acetic acid in each 
of the samples was determined by 
titration with 0.1 N standard sodium 
hydroxide solution with 1% phenol- 
phthalein as indicator. Twenty ml. 
of 95% ethanol were mixed with 
each ketone sample prior to titration 
to ensure the presence of only one 
phase (34), Each sample was corked 
as soon as it had been taken and was 
analyzed immediately following com- 
pletion of the run. 

All the experimental work was car- 
ried out at room temperature. The 
temperature of each of the phases 
both entering and leaving the column 
during each run and the room tem- 
perature during each run all ranged 
between 74.1°F. and 82.2°F., and 
averaged 77.5°F. 

Equilibrium data for the system 
methyl isobutyl] ketone-water-acetic 
acid were determined at 79°F. The 
results obtained agreed with those of 
Fleming and Johnson !9) for slightly 
lower temperatures. However, the 
effect of temperature is slight (2°). 

Caleulations. The average velocity 
of rise of the ketone drops was de- 
termined for many runs from the 
individual measurements of drop 
velocity. The ketone phase holdup 
was calculated from this average 
velocity by the method of Johnson 
and Bliss '*). The following equa- 
tion* was used: 


Li 
H = (100°; ) (1) 
u 

The water phase samples obtained 
by internal sampling showed no 
noticeable entrainment of the ketone 
phase. On the other hand, the ketone 
phase samples contained approxi- 
mately 50% by volume of the water 
phase when the glass sampling tube 
was used, and about 15% by volume 
of the water phase when the stainless 
steel sampling tube was used. The 
presence of a water phase in each 
ketone sample resulted in the trans- 
fer of acetic acid continuing after 
sampling took place. The volume of 
each phase of each ketone sample was 
measured before the respective con- 
centrations of acetic acid were deter- 
mined by titration. From these data, 
and from the known concentration of 
acetic acid at the sample point in the 
column obtained from analysis of the 
water phase sample, the original con- 
centration of the ketone drops could 
be calculated by a material balance 


°A table of nomenclature appears at the end 
of this paper. 
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represented by the following equa- 
tion (14): 


Vw (Cwt — Cwi) + Vice 
Cy = > (2) 
Ve 
When the equation was applied 


smoothed values of c,; were used. 
The trial and error procedure used by 
Gier and Hougen ‘!") in applying the 
material balance was not necessary in 
the present work. Any small changes 
in V,. and V, between the time of 
sampling and the time of analysis 
were neglected. 

The great majority of the ketone 
drops were observed to ascend the 
column continuously with a slight 
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zig-zag motion. (Gier and Hougen ‘!4) 
mention that their ether drops rose 
“erratically”.) The change in concen- 
tration of the ketone phase between 
any two levels in the column should 
provide a reasonably accurate measure 
of the mass transfer across the total 
interfacial area between those levels. 
The same statement can not be made 
with respect to the water phase be- 
cause of mixing taking place in that 
phase. The water phase concentration 
profiles resulted, not only from mass 
transfer through the interface, but 
also from such mixing. A hypothetical 
concentration profile for the water 
phase was calculated for each of 
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Figure 4—Ketone nozzle I. 
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Figure 5—Detailed drawing of ketone nozzle II. 


many of the runs. This profile would different equations, one based on the 
have been obtained under experi- change in concentration of the water 
mental conditions which produced phase between the two ends of the 
the measured ketone profile if the column, and the other based on the 
water phase concentrations were the corresponding change in_concentra- 
result only of mass transfer through tion of the ketone phase. These equa- 
the interface, or, in other words, if unions were 

true countercurrent flow 
the column. 
files were 


existed in a ee (3) 
These hypothetical pro- ; 7 
calculated from the 
smoothed ketone profiles by means of 
material balances written between the 


bottom of the 


Ni = LiS (cx, Ck?) (4) 


In general the value of N, was 
slightly different from that of N,. For 
use in further calculations N, and N, 
were averaged to give N. 


column and_ various 


heights above the nozzle tips. 


The rate of material transfer across 


the total interface in the column was As a partial check on the quality of 


calculated in |b. moles/hr. by two the experimental work the “percent 
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deviation” ‘3. +14) was calculated for 


each run: 


Bom Nw — Ne 
Percent deviation = «x 100% (5) 
Values of the overall capacity co- 
efficient, K,a, were calculated by 
means of the following equation ‘!4); 


2 


K,ah f dcx 
- = (6) 
Lx J (cx* — Ck) 
1 
The right hand side of the equation 
was evaluated graphically from the 
measured ketone phase and water 
phase concentration profiles. 
Smoothed values of the concentra- 
tions in the two phases were used. 
H.T.U. values were obtained from 
the mass transfer coefficients by 
dividing the latter by the respective 
flow rates of the ketone phase and in- 
verting the results (14 18); 


WE oe 7 
(ra. 3. tee oe (7) 
For purposes of comparison values of 
(H.T.U.),, were calculated on the 
basis of the assumption of true coun- 
tercurrent flow of the phases. Equa- 
tion 7.was applied to overall transfer 
co-efficients obtained by means of the 

following relationship ‘!*): 

K . 8) 

en Sh (cu* — ce), ; 
As will be indicated later in this 
paper the concentration end effect in 
the ketone phase at the continuous 
phase inlet to the column is repre- 
sented by the vertical distance AB in 
Figure 8. This end effect is believed 
to be due to an erratic jumping mo- 
tion of the ketone drops at the inter- 
face. It results in the part CD of the 
total concentration end effect in the 
water phase represented by distance 
CE in Figure 8. The remainder of this 
distance, DE, is believed to be the 
result of mixing, or circulation in the 
continuous phase by virtue of which 
the rising drops somehow carry with 
them continuous phase from near the 
bottom of the column and so cause an 
effective dilution of the continuous 
phase throughout the column for 
runs of the normal type previously 
defined. The amount of this dilution 
at the column top is represented by 


DE. 


6, the end effect due to motion of 
the ketone drops at the interface, was 
calculated as follows: 

B = LiS (concentration difference AB) 
(9) 
= LS (concentration difference CD) 


II 


y, the end effect in the continuous 
phase due to mixing, was obtained as 
follows: 


W =L,S (concentration difference DE) (10) 
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Equ: itions 9 and 10 as they have been TABLE 1 


given here refer to runs in which the EXPERIMENTAL CONDITIONS ; dl 
acetic acid was transferred from the spines ; : = r a 
water into the methyl isobuty]! ketone. Run Flow rates Temperature? Linear — Velocity, Holdup Column 
Where the reverse direction of ex- number! f/f. a . — - height? 
°"/ ° r. cetone in e ear . & 
traction was used § was zero and y : nozzle tips drop rise 
was calculated by Equation 10 with 
“concentration difference DE” re- Lw Lx u h 
placed by the total abrupt change in oF. — ft./hr. C; fe. 
acid concentration of the water phase ats : : : 
at the interface. 43 11.4 | 60.5 79 0.271 984 8.1 3.77 
; in 31 21.2 | 60.5 0.271 900 6.7 3.79 
Results. Forty-four runs were made 32 2121 60.5 74 0.271 3 78 
1 all. Some of these were discarded 33 20.1 | 89.1 76 0.266 865 10.3 3.80 
usually because the over-all material - = 7 = : a ; 275 1072 4.1 3 Poe 
a «ie . ial aaa ta S 7 7 eae aoa 
balance was poor. Some w ere repeats 39 42°7 | 18.3 72 0 286 1100 17 372 
done for various reasons including 23 42.5 | 35.8 0.280 1070 3.4 3 81 
checking the reproducibility of re- 18 42.5 | 37.5 79 0.293 1112 3.4 3.83 
sults and the effect of sampling on = 2 5 . 5 77 ¢ = = é 7 3 . 
Zo - ) 2 -44 3 
the steady state operation of the 29 42.5 | 605 0 271 937 65 3 81 
column. Tables 1 and 2 describe 17 27 42.5 | 89.1 75 0.266 810 11.0 3.80 
runs in which the normal direction of 28 42.5 | 89.1 = 0 _ . - 3.74 
extraction was used and five in which 30 85.0 | 60.5 a 0.271 930 0.9 3.79 
; 38 85.0 | 60.5 81 0.271 3.91 
the reverse was used. Table 1 lists the 41 | 84.9 | 89.2 78 0.266 754 11.9 3.76 
main experimental conditions under 
which the runs were carried out. bef as = ; 82 se - 
. Jc é. J. i | Zi 3 
ae Tes > rer-< ¢ =te i - - » 
Table 2 gives the over-all transfer 37 42.7 1| 605 80 0 271 3 Of 
data for each run. Tables 3 and 4 40 42.5 | 89.2 78 0.266 3.70 
show smoothed values of the point 42 84.9 , 60.5 80 0.271 3.80 
concentrations within the column for 
the water and ketone phases respec- 1In _ ae five runs only, the direction of extraction of acetic acid was from the ketone into the 
water phase. 
tiv ely. Fi igure 8 shows the concentra- 2Average temperature of aqueous phase, ketone phase, and room. 
tion profiles for Run 33, a normal 3Distance from nozzle tip to interface. 
run, and Figure 9, for Run 40, a re- Tin 2 
verse run. The dashed concentration OvER-ALL MATERIAL BALANCES 
curve in each figure represents the — =a 
water phase concentration profile ; : +. . ) 
ih eek a : ned Run Acid concentration Transfer rate Percent 
which would correspond to the auniher’ Deviation 
measured ketone profile for the con- lb. moles/ft.* X 105 Ib. moles/hr. X 10° 
dition of no circulation in the water : = 
phase. : 7 water ketone 
Average values of K,a and 
(H.T.U.),, were calculated by the : , 
° e ° ¢ 1 yu yu 
graphical integration procedure m = 7 N N N 
(Equations 6 and 7) applied to the | Cwr Cw2 Ck2 Cky a ~ — X 100 
whole column. Also obtained were ; 
¢ Vv > 4 ve s ~ 9 son ~ c ; a a - = e " 
average values of these quantities 438 50.40 12.92 515 12.18 5 24 0 38 
with the end effect in the ketone 31° 50.60 | 24.00 | 7.84 | 16.85 6.81 3.08 
phase eliminated. (This is the only 328 50.80 | 22.70 | 6.39 | 16.03 7.24 1.93 
end effect included in the average 33° 50.60 | 16.60 | 6.42 | 13.96 8.30 1.08 
ralues of K 1(H.T.U f a 13> 51.00 | 32.90 6.22 22.10 8.80 5.56 
values of Kya and (H.1.U.)o, for the 158 51.00 | 33.90 | 6.45 | 22.20 8.53 0.59 
whole column.) The end effect was 398 | 49.40 | 42.20 6.24 | 23.60 3.84 =e ae 
eliminated by stopping the graphical rea 3 00 7 50 oo a . 7.38 +.75 
: : ‘ b : ? 51.90 6.70 6.21 | 23.30 7.90 1.01 
integration at the last point below the ; : 
es | he 248 49.00 | 27.90 | 6.34 | 20.80 10.87 2.30 
interface at which samples of the 25> 44 40 | 26.90 6.13 | 18.00 9 04 4 54 
phases were withdrawn. Table 5 lists 29» 49.20 29.30 7.06 | 21.30 10.50 1.90 
the results of these calculations. In- 27» 52.30 | 23.30 6.31 | 20.30 15.22 “0.98 
: 284 53.30 | 22.50 6.08 | 20.40 15.89 2.39 
cluded also for purposes of compari- 302 51.50 | 39.20 . 36 | 24 es 12 82 0 30 
son are values of (H.T.U.),, based 388 51.60 | 39.10 | 7.59 | 25.10 13.03 —-0.23 
on the use of a logarithmic mean 41° | 49.90 | 30.90 | 5.51 | 23.90 19.95 1.50 
driving force calculated from the in- 
; let and outlet concentrations. 
Cw? Cw Cky Ck2 
Values of (H.T.U.),, were calcu- 
= 444! dé 5 aia 5 3 — 3 
lated for a few of the runs for short i: : 73 ; ba 3 as : o 7 on ; Mr 
lengths of the column by the graph- 374 0.17 | 11.02 | 12.40} 4.83 5.66 1.24 
ical integration procedure. The re- 40° 0.16 13.85 | 11.82 21 7.20 -1.53 
° da 7 ) ( on ? 
sults for Runs 30 and 35 are shown in 42 0.04 | 7.82 | 14.54) 3.49 8.16 1.22 
Table 6. 


1In the last five runs only, the direction of extraction of acetic acid was from the ketone into the 


The end effects § and y calculated water phase 
. ° " . : 2Internal sampling from nozzle tips to interface. 
7 € ‘ > r 
by Equations 9 and 10 are given in 3Internal sampling from interface to nozzle tips. 
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TABLE 3 
CONCENTRATIONS IN THE WATER PHASE! (SMOOTHED VALUES) 


Concentration of acid in water phase, lb. moles/ft.* X 10 


Run at distance from nozzle tips, ft 
number? 
0.1 0.5 1.5 25 3.5 h? 
43 13.0 13.1 15.1 21.2 29.3 30.6 
13.1 13.6 16.4 25.2 37.1 41.8 
31 24.8 26.6 28.6 33.6 37.9 38.7 
24.8 27.5 31.6 7.7 14.3 46.1 
33 7.2 19.2 21.5 242 32.7 34.0 
17.6 20.3 24.9 31.0 40.3 44.7 
39 42.5 13.6 45.7 47.2 48.1 48.3 
42.7 44.3 46.8 48.1 49.0 49.2 
23 38.0 40.1 43.9 47.0 49.0 49.2 
38.2 40.8 45.7 48.5 50.3 50.8 
24 28.3 30.5 35.4 40.3 44.4 45.6 
28.5 31.1 36.7 1.7 45.2 46.0 
27 24.3 27.6 33.8 39.2 43.8 45.1 
24.4 27.9 34.7 40.4 45.5 46.9 
30 39.7 41.6 45.0 47.4 49.6 50.1 
39.9 42.1 46.0 48.6 50.5 50.8 
41 32.0 35.5 41.2 44.8 47.7 48.3 
32.0 33.7 41.4 45.1 48.0 48.6 
44? 15.6 14.9 12.8 10.4 Lah 6.8 
15.0 13.1 9.5 6.0 1.9 0.3 
35 7.8 faz ae 4.1 2.6 2% 
7.4 6.0 3.6 1.9 0.7 0.2 
37 10.7 9.9 7.8 S.7 55 ae 
10.4 8.6 a7 3.3 i 0.2 
40 13.7 12.9 10.4 va 4.7 4.0 
13.4 12.0 8.3 4.6 0.9 0.2 
42 7.6 6.9 S41 5:3 1.4 0.9 
1.2 6.5 4.4 2.3 0.6 0.0 


' In each case the upper figure in the table is the value measured. 


[he lower figure was 


calculated by material balance on the basis of hypothetical true countercurrent flow of 


the phases. 


2 In the last five runs only, the direction of extraction of acetic acid was from the ketone 


into the water phase. 


3 Distance from nozzle tips equal to column height (Table 1). 
column are those just below the interface. 


The concentrations in this 


TABLE 4 
CONCENTRATIONS IN THE KETONE PHASE (SMOOTHED VALUES) 


Concentration of acid in ketone phase, Ib. moles/ft. X 10% at 


Run distance from nozzle tips, ft. 
number'! 

0.1 0.5 1.5 2.5 3.5 h?2 
43 5.3 5.4 5.9 von 9.7 10.7 
31 8.2 9.4 10.8 12.0 15.1 15.7 
33 6.6 7.2 8.3 10.0 11 8 12.6 
39 7.6 11.6 7.3 20.3 22.4 23.2 
23 8.2 3 16.6 20.1 22.3 i273 
24 6.8 8.6 237 16.3 18.7 19.1 
27 6.8 8.6 11.4 14.2 16.9 mt 
30 8.3 11.4 16.8 20.3 22.6 23.1 
41 6.7 10.0 15.4 18.9 2.7 an.8 
44! i. 11.4 9.9 8.7 7.0 6.5 
35 13.2 11.4 8.6 6.5 5.1 4.6 
37 12.0 10.8 8.5 6.9 5.4 +.9 
40 11.6 10.9 9.2 7.4 a2 xo 
42 14.0 12.6 9.7 734 4.8 4.2 


' In the last five runs only, the direction of extraction of acetic acid was from the ketone 


into the water phase 


? Distance from nozzle tips equal to column height (Table 1) 


column are those just below the interface 


Che concentrations in this 





Table No end effect due to agita- 
tion at the interface was found for 
the reverse runs 

Figure 10 is a photograph of a por- 
tion of the column under the operat- 


120 


ing conditions of Run 39. The ap- 
pearance of the drops in Figure 10 
was typical for normal runs. Most of 
the drops were ellipsoidal in shape 
and quite uniform in size. Judging 


from the known sizes of the sampling 
probes the horizontal axes of these 
drops averaged approximately 0.16 
in., and the vertical axes 0.09 in. for 
the conditions of Run 39. A number 
of small drops was always present 
with the larger drops in the column. 
These small drops were spherical in 
shape and varied from approximately 
0.01 in. to 0.05 in. in diameter. They 
were estimated to comprise about 
10% of the total number of drops, 
but would, of course, account for 
much less than this percentage of the 
interfacial area. Some of these small 
drops were probably removed with 
the water phase. Others, larger, rose 
with the ellipsoidal drops to the inter- 
face. Although no very accurate size 
measurements were made in runs 
other than 39, the size of the ellipsoi- 
dal drops was observed to remain 
approximately constant from run to 
run. No coalescence of drops was 
observed in normal runs. 


However, during reverse runs 
many of the drops coalesced within 
a short distance of the nozzle, and, 
before the interface was reached, 
large liquid masses with changing, 
irregular shapes resulted from con- 
tinued coalescence. S:milar behavior 
was noted by Johnson and Bliss ‘'%), 
A few small spheric: al drops were pre- 
sent in the reverse runs also. No 
coalescence of the small drops was 
observed. 


During normal runs the drops did 
not coalesce immediately on reaching 
the interface. From two thirds to all 
of the interface was covered with at 
least a single layer of drops, about 
one half of the interface usually was 
covered with at least two layers of 
drops, and about one third of the 
interface sometimes was covered with 
three layers of drops. The exact con- 
figuration of these layers depended on 
the water and ketone flow rates. The 
layers of drops were about half in the 
ketone phase and half in the water 
phase. At the middle of the interface 
considerable turbulent agitation of 
the drops was noticed. An accumula- 
tion of drops which had arrived from 
the column below would push up into 
the ketone phase above the interface 
and then subside by pushing drops 
radially toward the wall of the 6 in. 
inside diameter section of the column. 
As the drops moved outward they 
replaced drops that had coalesced. 
This spouting effect was most pro- 
nounced at high ketone flow rates 
and was almost absent at low ketone 
rates. Similar phenomena are des- 
cribed by Geankoplis and Hixson ‘"), 


During the reverse runs the liquid 
masses of ketone and the uncoalesced 
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drop: reached the interface and then 
coalesced immediately without any 
of the turbulent agitation noticed in 
the normal runs. 


Discussion 

Holdup. The ketone phase holdups 
listed in Table 1 were approximately 
directly proportional to the ketone 
phase flow rates. Similar results have 
previously reported, for ex- 
in Reference 14. 


been 
ample 


Contamination of the methyl iso- 
butyl ketone. The ketone phase re- 
maining at the end of Run 38 and 
containing most of the dissolved im- 
purities accumulated since the be- 
ginning of the work was concentrated 
by evaporation to a dark red residue 
which amounted to 3.4% of the 
original volume. Experiments showed 
that the equilibrium distribution of 
acetic acid between this concentrated 
residue and water was the same as 
that obtained earlier for acetic acid 
between technical grade methyl iso- 
butyl ketone and water. Furthermore, 
several values of K,a obtained by the 
logarithmic mean driving 
forces could be compared with those 
obtained by Johnson and Bliss (') 
under approximately similar operat- 
ing conditions. In most cases agree- 
ment was good. The presence of the 
impurities from the tubing does not 
appear to have affected the results 
reported in this paper. 


use ( of 


Sampling technique. The rate of 
sampling used, 10 ml./min., was ap- 
proximately 15% of the lowest 
aqucous phase flow rate of 11.4 
ft.*/hr.ft.2, and approximately 8% of 
_ next highest flow rate of 20.1 

‘/hr.ft.2. The corresponding figures 
co the ketone phase were 9% for 
pe fr.°/hr.ft.2 and 5% for 35.8 

/hr.ft.2. These percentages are all 
higher than those which applied dur- 
ing the ngeny done by previous 
workers (1. 2.3. 414), However, as can 
be seen from Table 1 only a few of 
the runs were made at these low flow 
rates. 

Samples of the phases leaving the 
column were obtained while internal 
sampling was being carried on. These 
samples were analyzed with the re- 
sults for a few of the runs shown in 
Table 7. The first sample for each 
run is believed to correspond to 
steady state operation of the column 
undisturbed by internal sampling. 
Internal sampling seems to have no 
significant effect on the concentra- 
tions of the phases leaving the column, 
even when the flow rate of one of the 
phases is low, as was the case for the 
Table 7. By con- 


runs in inference 
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Figure 6—Teflon caps and ‘tip patterns for ketone nozzle II. 


centrations throughout the column 


would also be unaffected. 


Purging or sampling through the 
sample probes was carried on almost 
continuously while each of the runs 
was in progress. Probably the main 
effect was to disturb the steady state 
by changing the flow rate of each 
phase at the location of the sample 
probes. Thus in Run 43 with the 
probes near the interface the flow 
rate of 11.4 ft.°/hr.ft.2 for the water 
phase would be 

11.4 — (0.15) (11.4) = 9.7 ft.8/hr.ft.? 
throughout most of the column. With 
the probes near the nozzle the flow 
rate would be 11.4 throughout most 
of the column. Thus the effect of 
sampling in this run was that the 
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profiles obtained were approximately 
those applying to an average con- 
tinuous phase flow rate of 10.5 

hr.ft.. where the flow rate was 
allowed to vary about 8% either way. 
The effect on the concentration pro- 
files was probably small, and this 
conclusion is supported by the data 
of Table 7 


It should be emphasized that Run 
+3 represents the very worst condi- 
tions of the sampling operation. In 
the reverse runs the sampling rate 
never exceeded about 7.5% of the 
flow rate of the respective phase and 
did not exceed about 4% except in 
In the normal runs it 
except in five in- 


two instances. 
did not exceed 5% 


stances. 
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DISPERSED PHASE SAMPLING TUBE 
INTO FLARED OUT 


PRESS FITTED 





1/8" TUBE 


SECTION. 


"0.D. x 0.020" WALL 


2" 
16 


CONTINUOUS PHASE SAMPLING TUBE. 


NOTE! BOTH TUBES CONSTRUCTED FROM TYPE 
304 S.S. SEAMLESS TUBING. 





Figure 7—Details of ketone and water sampling tubes. 


For normal Run 31 (L, = 21.2, 
L,, = 60.5) a sampling rate of approxi- 
mately 10 ml./min. represents respec- 
tively 8% and 3% of the flow rates of 
the phases. The internal samples were 
taken by moving the probes from 
sample point to sample point starting 
near the nozzle and moving toward 
the interface. In Run 32 the same 
experimental conditions were used 
except that the inlet ketone phase 
concentration was somewhat lower 
than in Run 31. The sampling was 
begun at the top of the column and 
successive samples were taken nearer 
spray nozzle. If 
sampling affected column operation 
appreciably, sampling in the different 
directions would probably 


and nearer to the 


produce 
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different concentration — profiles. In 
fact, although the concentration pro- 
files of Run 32 were displaced down- 
ward from those of Run 31 because 
of the lower ketone inlet concentra- 
tion, the profiles were of the same 
shape as in Run 31 throughout their 
lengths. Furthermore, both — runs 
showed roughly similar end effects 
(Table 5) of the two types discussed 
earlier. Thus it is concluded that the 
sampling rates were not excessive 
since the direction of sampling did 
not appreciably affect the concentra- 
tion profiles. A similar comparison at 
higher flow rates of the phases (Runs 
27 and 28) produced similar results. 

The sampling rates used and the 


presence in the column of the 


sampling tubes did not visibly affect 
the rise of the drops through the 
column except for a momentary hesi- 
tation in their upward course in the 
immediate vicinity of the ends of the 
probes. 

Samples of each phase were taken 
at various locations in the cross sec- 
tion of the column at a constant dis- 
tance from the nozzle tips. These 
samples were of constant concentra- 
tion for each of the respective phases 
except for deviations of one or two 
per cent, probably due to the analy- 
tical procedure used. 


Circulation in the continuous 
phase. The concentration of the con- 
tinuous phase was not constant over 
the column height in the present 
study. This result differs from that 
reported by Heertjes et al 5) but 
ee: with those of Geanpoklis et al 
(1, 2,8, 4), Therefore, at least in cer- 
we eee circulation is not 
so complete that the continuous 
phase is thoroughly mixed. 

The separation between the calcu- 
lated and measured water phase pro- 
files in Figures 8 and 9 is evidence of 
circulation in the continuous phase. 
In normal runs the calculated profile 
is above the measured profile because 
the water phase is diluted by the 
material circulated upward from the 
bottom of the column. However, in 
reverse runs the calculated profile is 
below the measured profile because 
relatively concentrated water phase 
is carried upward. Further quantita- 
tive information on the separation of 
the profiles is available in Table 3. 

The separation between the two 
water phase profiles depends in part 
on the amount and concentration of 
the material circulated, and on the 
flow rates of the phases. It might be 
expected that increasing the flow rate 
of the dispersed phase would have a 
marked effect on the separation be- 
tween the curves. However, com- 
parison of runs made at approxi- 
mately equal water flow rates but at 
different ketone flow rates showed 
that the amount of the separation was 
roughly the same. On the other hand, 
separation of the water phase profiles 
showed a marked decrease when the 
ketone rate was kept constant but the 
water rate increased: more water 
phase had to be diluted in the same 
time. The conclusion that the flow 
rate of the continuous phase is con- 
siderably more important than that of 
the dispersed phase in determining the 
amount of separation between the two 
continuous phase profiles is borne out 
by the results of both the normal and 
the reverse runs. 


Gier and Hougen ‘'#) discuss the 
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circulation in terms of a hypothetical 
flow pattern for the continuous 
phase. 

End effects. Evidence exists in 
some of the runs for an end effect at 
the entrance of the dispersed phase 
to the column. For example, see 
Figure 8 for Run 33. However, if 
present, this end effect is small and 
has been ~~. in this study. 
Kreager and Geankoplis describe 
a similar small effect. 

Inspection of Figure 8 and of 
similar plots for the other normal 
runs shows that for each of the 
phases abrupt changes in concentra- 
tion take place at = interface. In the 
reverse runs (e.g. I “igure 9), no dis- 
continuity appears in the dispersed 
phi ise profile at the interface, but one 
in the continuous phase profile per- 
sists. 

The abrupt change in concentra- 
tion of the ketone phase at the inter- 
face in normal runs is believed to be 
due to the agitation and spouting 
effect of the ketone drops prior to 
coalescence. The corresponding dis- 
continuity for the reverse runs is 
absent because the drops coalesce 
smoothly, immediately on reaching 
the interface. In Figure 8 the change 
of concentration of acetic acid in the 
ketone phase due to the agitation is 
shown as distance AB measured on 
the ordinate scale. This amount, 
multiplied by the flow rate of the 
ketone phase, gives the end effect , 
the rate of solute transfer into the 
drops during their stay at the inter- 
face. The acetic acid entering the 
agitated drops at the interface must 
be removed from the water. The cor- 
responding fall in the water phase 
concentration is represented by CD 
in Figure 8. By a material balance, CD 
times the flow rate of the water phase 
must equal the rate of solute transfer 
into the agitated drops, i.c. AB times 
the flow rate of the ketone phase. 
This is a statement of the relationship 
expressed in Equation 9. Point D can 
be located also by the material 
balance procedure used in establish- 
ing the water concentration profile 
shown as a dashed line. Values of { 
have been listed in Table 5 for many 
of the runs. 

The values of 6 include the effect 
of a variable driving force. This was 
eliminated by the following method. 
During the stay of the drops at the 
interface in the normal runs, mass 
transfer takes place as the result of 
a concentration difference: 


CF? = 


An average value of this driving force 
was obtained as follows. A value of 
the expression was calculated from 


o 
oO 


a 
oO 


BS 
°o 
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ACID CONCENTRATION , LB. MOLES/CU.FT. x 10° 
bad 
° 


BTM. 1.0 





Ly, =89.1 CU. FT./(HR.)(SQ.FT.) 
Ly =20.0 CU. FT./(HR.)(SQ.FT.) 


CALCULATED WATER 
PHASE PROFILE, .- 


WATER PHASE 
PROFILE 


INTERFACE — 


OwA 
O--B 


KETONE PHASE 
PROFILE 


2.0 3.0 4.0 


COLUMN HEIGHT, FT. 


Figure 8—Conceentration profiles for Run 33. 


the entering water phase concentra- 
tion and the leaving ketone phase con- 
centration at the top of the column. A 
second value was obtained from the 
concentrations measured just below 
the interface during internal sampling. 
(Smoothed concentrations were not 
used here.) The arithmetic average of 
the two values was calculated. This 
average driving force was then di- 
vided into the values of 6 previously 
obtained to give quotients: K,A, 
listed in Table 5. The result of plor- 
ting K,A against the ketone phase 
flow rate with the water phase flow 
rate as a parameter is shown in Figure 
11. Although the points scatter rather 
badly, relationships are evident be- 
tween K,A and both the ketone phase 
flow rate and the water phase flow 
rate. As would be expected from the 
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description given previously of the 
behavior of the drops at the interface, 
K,A increases with the ketone phase 
flow rate. K,A increases also with 
the water phase flow rate. This last 
result indicates that considerable dis- 
turbance exists at the top of the 
column as a result of the entry of the 
continuous phase even though the 
design of Blanding and Elgin “% was 
used for this part of the apparatus. 
As would be expected, and as Figure 
11 indicates the amount of this dis- 
turbance is higher at higher 
flow rates. 

Phe quantity of K,A is similar to an 
over-all capacity coefficient Kya. 
However, the area factor in K,A 
is the total interfacial area of the 
drops at the interface, plus that part 


water 


123 








TABLE 5 
AVERAGE V ALUES OF vay AND oF (H.T.1 


.)ok AND END EFFECTS 





\verage values of Kya and of (H.T.U. )ox End effects 
Flow rates Mass transfer 
coefficient for 
Results of graphical integrat Results Ib.moles/hr. X 10° turbulence 
Run ft.3/hr.ft.? using at 
number? log interface 
Whole column lo point bek iterface mean Turbulence Circulation 
sa I J 
Kya HT. Kya? H HT.) B y KyA 
ft ft ft.3/hr. 
43 11.4 60.5 36.6 1.65 30.7 1.97 2.70 117.4 588.0 0.126 
31 24.2 60.5 34.0 1.78 30.7 1.97 2.66 100.4 209.0 0.146 
32 212 60.5 30.5 1.98 24.9 2.43 2.87 158.3 223.5 0.213 
33 20.1 89.1 54.7 1.63 50.1 1.78 2.66 138.0 267 .0 0.162 
13 40.7 44.0 28.8 1.53 26.6 1.65 1.65 59.5 60.0 0.111 
15 40.7 44.0 27.0 1.63 23.9 1.84 ‘71 59.5 60.0 0.135 
39 42.7 18.3 14.1 L320 13.0 1.41 1.48 1 a $2.5 0.055 
23 42.5 35.8 25.4 1.41 2a 1.54 1.63 44.0 107.8 0.129 
18 42.5 37.5 25.9 1.45 22.6 1.66 1.66 59.7 122.8 0.171 
24 42.5 60.5 44.1 Lal 38.8 1.56 1.60 112.5 57.5 0.235 
25 42.5 60.5 38.5 Pan 33.6 1. 80 1.89 92.8 130.4 0.220 
29 42.5 60.5 42.3 1.43 35.8 69 1.58 119.0 37.6 0.256 
27 42.5 89.1 59.8 1.49 §1 2 74 1.70 274.0 97.1 0.427 
28 42.5 89.1 63.6 1.40 56.4 1.58 1.61 197.0 179.0 0.313 
30 85.0 60.5 48.4 1.25 40.3 1.50 1.36 100.4 66.9 0.375 
38 85.0 60.5 54.5 2.11 46.5 1.30 1.27 52.0 156.0 0.289 
41 84.9 89.2 86.5 1.03 70.2 27 1.06 197.0 1155 0.758 
44! 22.8 60.5 23.1 2.41 25.1 2.41 3.65 183.0 
35 42.8 35.6 ‘.2 2.353 14.2 2-3) 3.01 102.5 
37 42.7 60.5 24.3 2.49 24.3 2.49 3.01 135.0 
40 42.5 89.2 38.6 2300 38.6 2.31 2.95 201.0 
42 84.9 60.5 26.0 rise 26.1 2.32 2.24 93.9 
‘In the last five runs only, the direction of extraction of acetic acid was from the ketone into the water phase. No ketone end effect 
was observed. Hence for these runs only, Column 4 in the table is the same as Column 6 and Column 5 is the same as Column 7 
? Units: Ib.moles/(hr. )(ft.*)(Ib.moles per ft 
TABLE 6 [ABLE 7 
(H.T.U. ox For SHORT LENGTHS 01 \NALYSIS OF EXIT SAMPLES 
COLUMN 
Run Flow rates Time after initiation of Concentration of 
i allele number internal sampling! to acetic acid in 
Run seats tips ft (H.T.I ft.3/hr.f taking of exit samples phase leaving column 
mber : 
— lb.moles/ ft.* 
From lo ft wa éwcteiat 
308 0 0. 240 106 Lw | min. water ketone 
0.240 0.515 1.40 - x aa a seas 
0.515 0 865 1 a i) 11.4 60.5 8 0.01330 0 01220 
1 g 1 16 1 5() 50 0.01272 0 01218 
1.7 2.03 1.47 - Pa 
03 31 1.38 1 21.2 60.5 2 0.0244 0.0173 
> 31 9 62 1 35 33 0 0229 0.0168 
2 62 3.00 1 4) 53 0.0248 0.0166 
3.00 3.747 1.92 
» 47 2 mo 9 42.7 18.3 4 0.0422 0 0234 
786° 0. O7¢ 
nae er 26 0.0424 0.023 
57 0.0421 0 0236 
Arithmetic average first 10 values 1.46 ; 
Including purging of the sampling lines. 
35% 0 0.15 1.83 
0.15 0. 64 eae 
0.64 Pe » 46 
1.35 ; 2.28 2.52 . . . ; ; 
» 28 2 Od 29 of the area of the interface itself not phase. As indicated previously, point 
covered by drops. On the other hand, D may be located by balancing two 
Aeteh mentio averace 2 34 1, IN Kya, is the interfacial area per transfers: the acetic acid gained by 


® Extraction of 


ketone 


® Extraction of 


water. 


aceli¢ 


aCelle 


Point just below interface 


4 [nterface 
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top of column 


acid trom water 


icid from ketone 


to 


cubic foot of column. 

In Figure 8, point D would repre- 
sent the concentration of 
in the water phase just below the 
when the measured ketone 
phase profile existed, and for condi- 
of no circulation in the water 


acetic acid 
interface 


ohlons 


the ketone phase in moving between 
the bottom of the column and the 
point just below the interface, and an 
equal loss of acetic acid by a water 
phase flowing in true countercurrent 
manner from the point just below the 


interface to the bottom of the 
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Figure 9—Concentration profiles for Run 40. 


column and leaving at the measured 
exit concentration. Absence of circu- 
lation in the water phase would result 
in point D coinciding with point E 
on the measured water phase concen- 
tration profile. However, as the ke- 
tone drops rise, they carry upwards 
“atmospheres” (3) or “envelopes” ‘*) 
of the continuous phase from regions 
of relatively low acetic acid concen- 
tration. When the drops coalesce they 
discharge these “atmospheres” which 
tend to dilute the water solution of 
acetic acid at the interface and pro- 
duce the separation between points 
D and E. The product of DE and the 
rate of flow of the water phase repre- 
sents the number of moles per hour 
apparently transferred out of the 
water phase and into the ketone phase 
because of the circulation § effect 
(Equation 10). This quantity has been 


called y in this paper. 

Values of y are shown plotted in 
Figures 12 and 13 for the normal 
runs. In Figure 12 the dashed line is 
the best straight line through the 
points. The solid line is the best 
straight line through the points and 
the origin. Circulation in the con- 
tinuous phase increases only slowly 
with the ketone phase flow rate. 

A few data estimated from the 
work of Gier and Hougen \!") have 
been included in Figures 12 and 13 
for the purposes of comparison. These 
data were obtained with a column 
having a diameter to height ratio 
approximately three times that used 
in the present research. From Figure 
12, decreasing the diameter to height 
ratio seems also to decrease the circu- 
lation as suggested in Reference 14. 
However, the measurements of Gier 
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Figure 10—Drops in column, condi- 
tions of Run 39. 


and Hougen were made with a differ- 
ent system from that used in the 
present work so that further measure- 
ments are contemplated in which the 
same system is used in towers with 
various diameter to height ratios. The 
effect of length to diameter ratio on 
back-mixing in the gas phase of fluid- 
ized beds, an analogous problem, has 
been studied by Gilliland and Mason 
21, 22) 

Figure 13 shows that circulation is 
greatly increased when low water 
rates are used. The end effect due to 
circulation is large enough to produce 
a very large value of y although the 
low value of the water rate is used 
as a multiplying factor with this end 
effect in calculating y. 

A possible explanation of the results 
shown in Figures 12 and 13 may be 
obtained if it is postulated that a small 
change in the velocity of the drops 
to the continuous phase will produce 
a large change in the circulation. At 
low water phase flow rates, the hold- 
up of dispersed phase should be a 
little less than it is at high water 
phase rates for the same ketone flow 
rate. The resulting slightly lower de- 
interference 
drops, and consequent higher velocity 


gree of between the 


of the drops relative to the water 
phase, possibly may account for the 
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high values of y at low water rates 
in Figure 13. On the other hand, 
increasing the holdup of drops by 
increasing the ketone phase flow rate 
at constant water phase flow rate will 
result in increased interference be 
tween drops so that the velocity of 
the drops relative to the wate! phase 
is decreased [hus, as Figure 12 
shows, Increasing the ketone rate al 
constant water rate does not produce 
large increases im y 

The data for the series of runs 
numbered 39 to 27 in Lable | provide 


experimental evidence for the de 


126 


creased velocity of the drops relative 
to the continuous phase resulting 


trom the ketone rate being increased 


g 
it constant water rate. These runs 
were all made at a water rate of 42.5 
to 42.7 ft./hr. The table shows that 
the velocity of rise of the drops rela 
tive to the equipment decreased from 
1100 ftt/hr. at a ketone rate of 18.3 
tt.°/hr.tt.2 to 810 ft./hr. at a ketone 
rate of 89.1 ft.4/hr.ft.*. 

In the reverse runs the whole end 
effect in the concentration profile for 
the continuous phase is believed to be 


due to circulation. The end effect in 


the ketone phase profile was absent 
because the drops coalesced smoothly 
as soon as they reached the interface, 


Hence, to calculate y for the reverse 
runs the abrupt change in concentra- 
tion of the water phase at the inter- 
face was simply multiplied by the 
How rate of the water phase 

Figure 14 shows the dependence of 
w on the flow rates of the two phases 
for the reverse runs. The general 
pattern of the variations in yw are 
those already described for the normal 
runs. 


Capacity coefficients and 
(H.T.U.),,. values. The capacity 
coefficients and H.T.U. values shown 
in Table 5 agree reasonably well for 
runs made under the same flow con- 
ditions of the phases. An exception 
to this statement is found in the case 
of the pair of runs numbered 31 and 
32 previously referred to in the dis- 
cussion of sampling technique. Here 
the results for the whole column 
check fairly well, but the results ob- 
tained when the end effect was 
eliminated do not. 

For the normal runs values of K,a 
ire reduced and values of (H.T.U.),,, 
ire increased when the graphical inte- 
gration is stopped just below the 
interface instead of being continued 
to include the concentration of the 
ketone leaving the column and that 
of the water entering. This result is 
obtained because agitation at the 
interface accomplishes a considerable 
mass transfer in a short distance. For 
the reverse runs the concentration of 
the ketone drops just below the inter- 
face is the same as that of the coa- 
lesced ketone phase above the inter- 
face. Hence crossing the interface 
results in zero contribution to the 
irea under the curve used for the 
graphical integration. (c, is the 
ibscissa of this curve. Thus rhe 
results for the whole column are the 
same as those to the point just below 
the interface. In other words, there is 
no ketone phase end effect, and, 
therefore, no correction to apply to 
the mass transfer coefficients and 
H.T.U.’s. 

The results for (H.T.U.),, ob- 
tained by graphical integration over 
the whole column often differ 
markedly from the (H.T.U.),,, values 
calculated from a logarithmic mean 
driving force based on the inlet and 
outlet concentrations of the phases. 
Chis result is particularly evident in 
the normal runs at low water rates at 
the beginning of Table 5 and in all 
the reverse runs at the end of the 
table. For both sorts of runs circula- 

on has produced a marked separa- 
tion between the experimental and 
the calculated continuous phase pro- 
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those for the normal runs. . | 
: ‘ > 
Table 6 indicates that the value of 
(H.T.U.),,, for normal runs probably ol / 
remains fairly constant up the column. oO 20 40 60 80 i00 
ar } > “Tak trek we y 2¢ 
However, in the reverse Run 35 the WATER PHASE FLOW RATE, CU.FT./(HR.)(SQ.FT.) 
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- ee = tg 
interface, probably mainly as a result Figure 13—Variation of y with water flow rate. 
of decreased interfacial area caused 
by the drops coalescing. In other of 
the reverse runs the (H.T.U.),, in 
creased toward the interface in the 300 
bottom parts of the column, but 
apparently fell off again as the inter- © 
*e r¢ approached. The 1 “rns ° 
face was ap} roached The a a] oO W VS Ly AT CONSTANT L, OF 
sd ‘ = . pres- 
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ent work are reasonably adequate for 200 
studving circulation and end effects, ow 
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end effect due to the turbu- 
lent agitation of drops at the 
interface, as defined by 
Equation 9, lb.moles/hr. 


end effect due to circulation 
in the continuous phase, as 
defined by Equation 10, 
Ib.moles/hr. 


c integral between elevation 1 
and elevation 2 in the col- 

J umn. 

1 
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Waste Treatment at Sun Oil’s 
Sarnia Refinery 


HAROLD F. ELKIN? and GORDON R. HENDERSON? 


The design of the waste water treatment plant at the Sarnia refinery 10,000 bbl. cone roof tanks are avail- 
of the Sun Oil Co. Limited is described. Operating results for the sulphide able for ship ballast water disposal. 
stripping tower and the biological oxidation unit are given. Over a two The process sewer system carries 
month test period 91.0% of the sulphides and 95.5% of the phenol were : 


- the main oxygen-consuming contami- 
removed and the average cost of phenol removed was 0.45 dollars/lb. 


nents in the refinery, 


HE combined efforts of industry 
‘aa the International Joint Com- 
mission On The Pollution Of 
Boundary Waters to preserve the 
water quality of the St. Clair River 
below the Sarnia, Ontario area have 
received deservedly wide attention (!), 
That these efforts are meeting with 
success is exemplified by the con- 
tinued attempts to adhere to I. J. C. 
objectives by existing and new indus- 
tries. ; 

Design of Sun Oil Company’s Sar- 
nia refinery in 1952 presented a 
unique opportunity to construct pol- 
lution abatement facilities as an in- 
tegral part of the plant and to pro- 
vide for expansion of waste treatment 
ex quipment as required by operation 
experience. This paper will describe 
the development of this installation 
and discuss operating results obtained 
to date. 


Sewer and Separator Facilities 
Sun’s Sarnia refinery, which was 
designed to process 15,000 bbls./day 
of crude oil, contains completely in- 
tegrated facilities for control and 
treatment of plant waste waters. The 
original installation which was placed 
in operation in 1953 included an ex- 
tensive segregated sewer and separa- 


tor system which served as a basis for 
current American Petroleum Institute 
recommended standards for new re- 
fineries ‘?). Independent drainage sys- 
tems were installed for the process 
waste waters, contaminated storm 
drainage, spent cooling water, uncon- 
taminated storm drainage and sanitary 
sewage. A sixth sewer system for 
sulphide-bearing process waters was 
added in 1955. The process wastes, 
contaminated storm waters and spent 
cooling waters are each provided with 
individual oil-water separators of API 
design for removal of oil and settle- 
able solids before discharge to the 
efHuent sewer. The total flow of the 
combined waters in the final effluent 
sewer averages 14 million gallons per 
day. 

The contaminated storm water 
separator unit consists of two rein- 
forced concrete parallel single-section 
chambers, each 20 ft. wide x 160 ft. 
long x 10 ft. deep. The cooling water 
separators consist of two parallel 
single-section chambers, each 20 ft. 
wide x 180 ft. long x 10 ft. deep. Two 
1000 bbl. cone roof tanks are provided 
for slop oil recovery and three 70 ft. 
wide x 100 ft. long x 6 ft. deep earth- 
wall storage basins are available for 
separator sludge dewatering. Two 


iManuscript received April 5, 1957. 
“Sun Oil Co. Limited, Sarnia, Ont. 


3Manager, Sarnia Refinery, Sun Oil Co. Limited, Sarnia, Ont. 


Contribution from the Sun Oil Co. Limited, Sarnia, 


Ont. 


Presented at the Third Ontario Industrial Waste C onfe rence, June 11-13, 


1956, Honey Harbour, Ont. 
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averaging 200 
GPM or somewhat less than 300,000 
gallons per day. This flow contains 
over 95% of the phenols, sulphides, 
mercaptans, emulsified oils and other 
water-soluble pollutants. These waters 
come from contact with oil in the 
plant and include process condensate 
drains from vessels and towers, pump 
gland leakage and desalting water. A 
10 ft. wide x 40 ft. long x 8 ft. deep 
concrete single-section oil-water sepa- 
rator was provided for primary treat- 
ment of this process waste water. It 
was anticipated that these process 
waters would require secondary treat- 
ment before discharge, and their 
segregation was considered a practical 
initial step. The general flow diagram 
for segregation of refinery waste 
waters, from which the Sarnia sewer 
and separator system was developed, 
is shown in Figure 1. An over-all view 
of the Sarnia refinery separator area 
is shown in Figure 2. The effective- 
ness of segregation of wastes at their 
source in aiding subsequent oil-water 
separation has been established in pre- 
vious investigations, ‘3) and the above- 
described Sarnia refinery drainage 
system has continuously produced 
final net effuent with an oil content 
of less than 5 ppm. 

Secondary treatment for removal 
of soluble materials is accomplished 
by flue gas stripping for removal of 
sulphides followed by biological oxi- 
dation in an activated sludge unit for 
reduction of phenols and other oxy- 
gen-consuming contaminants. 
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2—An over-all view of the Sarnia refinery separator area. 
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Figure 3—Flow diagram for sulphide stripping tower—1955. 


Sulphide Stripping Tower 

The flow diagram for the sulphide 
— tower installation is shown 
in Figure 3. This unit utilizes waste 
flue gas to remove objectionable sul- 
phides and mercaptans from the pro- 
cess water system and to add temper- 
ature to the feed water to the subse- 
quent biological treatment unit. Ap- 
proximately 75 GPM of the process 
waste waters are collected in an 
underground concrete oil separating- 
surge sump and are then charged to 
the top of the 5 ft.-6 in. diameter by 
32 ft. high stripping tower through 
six 1} inch full-cone spray nozzles. 
The injected feed water flows down- 
ward against a continuous rising 
column — of approximately 75,000 
SCFH of 950°F. flue gas from the 
refinery catalytic regenerator to re- 
move over 90% of the sulphides from 
the feed water. Operation data for the 
stripping tower during recent months 
is shown in Figure 4 and summarized 
in Table 1 


This stripping tower is unique in 
refinery practice in several _ 
Because flue gas was available only a 
the low pressure of 2 psig., the entire 
tower and gas-piping system had to 
be designed for a maximum pressure 
drop of 1.5 psig. A tower design was 
therefore developed consisting only 
of a distribution header with the six 
spray nozzles and a demisting section. 
No tower packing or trays were 
installed. The internal walls of the 
tower shell were lined with 2 inches 
of gunnited high alumina-vermiculite 
cement for corrosion resistance. Be- 
cause of the intimate contact of the 
sprayed feed water with the flue gas, 
sulphide removal has been satisfactory 
even at low temperatures as noted in 
the data of Figure 4. It has therefore 
not been necessary to augment the 
Hue gas with stripping steam to main- 
tain performance. 

The stripping tower desulphurized 
bottoms are heat exchanged with the 
sulphide-rich charge water and then 
flow with the other process waste 
waters to the process separator for 
further oil removal before additional 
treatment. 


Biological Oxidation Facilities 

Sun Oil Company initiated opera- 
tion of commercial-scale biological 
oxidation for refinery waste waters at 
the Sarnia plant in 1953. A modified 
activated sludge unit was installed to 
oxidize phenols and other organic 
contaminants in the process waste 
waters. The activated sludge process 
was selected over competitive 
methods of treatment because of 
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relat've simplicity of construction, TABLE 1 
low investment cost, minimum space SULPHIDE STRIPPING TOWER OPERATION DATA-SUMMARY* 


yon srs SARN SFINERY — Sun O1L Co. LIMITED 
requirements, availability of low-cost Sarnia REFINERY — Sun On Co. Liarrer 


electric power and anticipated high 

ope ation efficiency. The utilization Stripper Influent Stripper Effluent Per 
of sludge- recirculation in the aeration > he een eecl, . —_ ' 
basin for metabolic oxidation (endo- Nien Min. | Avg Max Min Av, AVE 
genous respiration) of excess bac- - - ees 
terial growths, thereby eliminating Flow — GPM : 
separate sludge digestion equipment, 
represented another advantage of the 
activated sludge process. Sulphides — ppm 1000 424 705 . 91.0% 


Temperature — °F 109 68 


When Sun Oil developed plans for 
the biological treatment unit, helpful *Data based on March-May 1956 operation. 
information was made available by 
the Dow Chemical Company ‘*) and 
Imperial Oil Limited \.®) who had 
plant and pilot-scale experience re- 
spectively with biological removal of 
phenols. Construction of the bio- 
oxidation unit was expedited by the 
use of the process sewer and separator 
system, in which the phenolic waters 
had already been segregated, and by 
utilizing other existing equipment 
and piping already in place. As shown 
on Figure 5, the aforementioned sul- 
phide stripping tower and process 
separator are necessary initial steps in 
pons water conditioning. The 200 
GPM flow from the process separator 
is transferred by a self-priming centri- 
fugal pump, rated at 500 GPM and 35 
psig. discharge pressure, to an earth- 
wall oxidation basin. The aeration 
basin has top dimensions of 90 
ft. long x 65 ft. wide x 10 ft. deep. 
The bottom half of the side-walls of 
the basin are surfaced with a 2” layer 
of cement grout to minimize erosion. ibis 
The fresh feed water and recycle Trt ea 
(mixed liquor volume totals 500 
GPM) enter the bottom of the basin 
through a nest of 40 mixing eductors. 
A rotary blower, rated at 500 cfm at 
psig., supplies air to the eductors 
providing oxygen and intimate mix- 
ing. The overflow from the aeration 
basin flows through a 14” peripheral 
slit-pipe flume to one section of the [enor =— 
contaminated storm water separator 
originally designed .for oil removal aika 
service. This separator is equipped ilies Fama 
| | 
| 


SULFIDES 
PPM 





TEMPERATURE 





Figure 4—Sulphide stripping tower operation. 
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serves as a secondary settling tank for nae, 
bacterial sludge. A positive displace- soto" x45" 6" aera 
ment sludge pump, rated at 125 GPM 
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terial floc to the oxidation basin. 2 Care 
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finery discharge sewer to the St. Clair 
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Figure 6—View of the oxidation basin during construction. 





Figure 7—View of the aeration basin in operation. 


TABLE 2 
BIOLOGICAL OXIDATION PLANT OPERATION DaTA-SUMMARY* 
SARNIA REFINERY Sun O1 Co. LIMITED 


Plant Influent Plant Effluent Per 
Z Cent 
Removal 

Max Min Avg Max Min Avg Avg 
Flow GPM 250 150 200 250 150 200 
pH 8.9 8.3 8.8 8.2 a 7.8 
Sulphides — ppm 76 0 22 0 0 0 

Sulphides — Ibs/day** 182 0 55 0 0 0 100% 
Phenol — ppm 56 14 41 4.8 0.4 ue 

Phenol — lbs/day** 135 34 98 12 1 4.5 | 95.5% 


*Data based « 


**Based on ave 
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yn March-May, 1956 operation. 
rage flow rate. 


Figure 6 shows a view of the oxida- 
tion basin during construction. Views 
of the aeration basin in operation are 
shown in Figure 7. 

The modified activated — sludge 
facilities were placed in operation in 
November, 1953, in time to treat the 
first waste waters from the refinery. 
Two tank truck loads of concentrated 
bacterial sludge from the Dow Chemi- 

cal Company waste treatment plant 
at Midland, Michigan, were obtained 
for initial seeding. Bacterial oxidation 
of phenol began within the first few 
hours. The approximate first cost of 
the new facilities was $40,000. This 
does not include the value of the pre- 
viously installed equipment including 
the process separator, collecting sys- 
tem and the secondary settling tank 
nor does it include the recently in- 
stalled sulphide stripping tower. 


Experience through the first two 
years of operation indicates that an 
average phenol content in the process 
waste feed of approximately 100 
lbs./day can be oxidized satisfactorily. 
Daily loadings as high as 200 lbs. of 
phenol have been removed with no 
loss in efficiency. Numerous mechani- 
cal difficulties and variations in feed 
water quality interfered with con- 
tinuous efficient performance of the 
big-oxidation unit during the first 
year of operation. High alkalinity and 
sulphide concentrations in the feed 
water have been corrected by the 
above mentioned flue gas stripping 
tower. Excessive sulphides formerly 
exerted a high chemical oxygen de- 
mand on the system and proved toxic 
to the bacterial colony. Efficient 
operating conditions include a_bac- 
terial sludge volume index of 4-6%, 
a temperature range of 70-90°F. and 
a pH range of 6.0-9.0. Experience in- 
dicates that the process will operate 
with a temperature as low as 60°F. 
but with reduced efficiency. Dissolved 
and emulsified oil concentrations as 
high as 150 ppm., phenol concentra- 
tions up to 200 ppm. and sulphide 
concentrations up to 50 ppm. have 
had no pronounced harmful effect on 
biological activity, although higher 
values can seriously injure the bac- 
terial sludge. High oil content in the 
feed water, although apparently not 
toxic to the organisms, is adsorbed on 
the biological floc and floats the mix- 
ture to the surface of the basin, there- 
by removing organisms from the pro- 
cess. It is believed that some oil is par- 
tially oxidized. Steam is introduced 
for temperature control purposes 
through 1” steel pipe coils with sub- 
merged exhausts into the oxidation 
basin. A_ direct-contact 2,500,000 
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BTt ‘HR water heater has been in- 
stalicd above the surface of the basin 
as a:) auxiliary source of heat. 


Ir has been found important to 
avo.d sudden or shock loadings in in- 
fluent characteristics. It has been ob- 
served that organisms can remain in- 
active for several days under adverse 
conditions, then regain their former 
activity with the return of normal 
teniperature and food supply. The 
hardiness of the bacterial colony is 
best exemplified by the fact that re- 
seeding of the unit from a neighbor- 
ing sewage plant has been required 
only once since original start-up, fol- 
lowing a scheduled refinery main- 
tenance shut-down. During other 
periods when refinery units were not 
operating and there was no con- 
tinuous flow of waste process water, 
mixed liquor was recirculated at the 
oxidation basin and chemical-grade 
phenol was used to continue the food 
supply. 

Air requirements of the aeration 
basin are supplied by 40 jet eductors 
mounted above the floor. Each aera- 
tor supplies 50 lb. of oxygen or a total 
of 2,000 Ib./day to the waste water 
to oxidize the average 130 lb./day 
(maximum 350 lb./day) of phenol, 
considerably in excess of the theo- 
rectical stoichiometric ratio of 2.4 to 
1.0. The aeration tank loading is 
equivalent to 4 lb. of phenol per 1,000 
cu. ft. of tank volume, and the applied 
air rate is equivalent to 0.1 cu. ft. per 
minute per sq. ft. of aeration surface. 
Total power requirements of the 
pumps and blower amount to 60 HP. 


When the above-mentioned pro- 
blems of temperature and_ influent 
characteristics were corrected, opera- 
tion of the biological treatment unit 
produced a continuously satisfactory 
efHuent. Recent operating data are 
presented graphically on Figure 8. A 
summary of average, maximum and 
minimum phenol contents of biologi- 
cal treatment unit influent, effluent 
and final combined refinery discharge 
are tabulated in Table 2. As evi- 
good operation is 
over 95% removal 


denced herein, 
characterized by 
of phenols. 

One of the remaining operating 
difficulties still to be corrected is the 
continuing use of the rectangular oil- 
water separator as the secondary bio- 
logical sludge settling tank. This 
equipment is over-sized and not of 
conventional design for bacterial floc 
recovery. Although oxidation effi- 
ciency is not seriously affected, a 
suitably- -designed replacement tank is 
under consideration. 
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Figure 8—Biological oxidation plant operation. 


TABLE 3 


ANNUAL Cost OF 
BIOLOGICAL OXIDATION 


Operations. <0 6s. cea 
Maintenance... 
Electric Power. 


Oe ae 


Total Annual Cost ‘ 
Cost Per Pound Phenol Removed... . 
Initial Capital Investment 


Annual operating and maintenance 
expense for the biological oxidation 
unit during the past year are sum- 
marized in Table 3. The unit cost of 
$0.45 /lb. of phenol removed is con- 
sidered moderate when compared 
with alternate chemical treating or 
incineration procedures. 

The two years of operation of the 
Sarnia unit indicates that the acti- 
vated sludge process is applicable to 
refinery waste waters. Despite early 
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PHENOL 
SARNIA REFINERY, SUN O1L Co. 


REMOVAL 1955 


LIMITED 


, 000 
, 800 
4, 200 
6,000 


Uw Ww 


$19,000 
$ 0.45 
$40,000 


operating difficulties, the equipmert 
has displaved good reliability when 
supplied with a consistently accept- 
able feed water and has provided a 
satisfactory solution to the plant 
organic pollution problem. 


Conclusion 

The opportunity of designing pol- 
lution control facilities as an integral 
part of plant operation equipment 
represents an effective approach to 
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treatment requirements. Installing the 
segregated sewer and separator system 
during initial Sarnia refinery con- 
struction made possible the subsequent 
rapid and economical installation of 
the biological oxidation unit to treat 
a limited volume of concentrated 
waste waters. 
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